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ABSTRACT 

We present the first results from our spectroscopic survey^ of the environments of strong gravi- 
tational lenses. The lens galaxy belongs to a poor group of galaxies in six of the eight systems in 
our sample. We discover three new groups associated with the lens galaxies of BRI 0952—0115 (five 
members), MG 1654+1346 (seven members), and B2114+022 (five members). We more than double 
the number of members for another three previously known groups around the lenses MG 0751+2716 
(13 total members), PG 1115+080 (13 total members), and B1422+231 (16 total members). We 
determine the kinematics of the six groups, including their mean velocities, velocity dispersions, and 
projected spatial centroids. For the newly discovered groups, we quantify these properties for the first 
time. For the other three groups, the increased membership allows us to make more robust estimates 
of the kinematic properties of the groups than previously possible. The velocity dispersions of the 
groups range from llOj^g™ to 470j^gQ" km s^^. The higher velocity dispersions (for the richer groups 
MG0751, PG1115, and B1422) are consistent with those of nearby X-ray luminous groups, while the 
others (for the poorer groups BRI0952, MG1654, and B2114) are more typical of nearby dynamically 
younger groups. The lens galaxy is the brightest member in fewer than half of the groups. In general, 
the brightest group galaxy is an early-type galaxy that lies off the center of the potential and occupies 
an orbit indistinguishable from the other group members. In at least three of the lenses — MG0751, 
PG1115, and B1422 — the group environment significantly affects the lens potential. These lenses 
happen to be the quadruply-imaged ones in our sample, which suggests a connection between image 
configuration and environment. Finally, our survey allows us to assess for the first time whether mass 
structures along the line of sight are important for lensing. We first show that, in principle, the lens 
potential may be affected by line-of-sight structures over a wide range of spatial and redshift offsets 
from the lens. We then quantify real line-of-sight effects using our survey and find that at least four 
of the eight lens fields have substantial interloping structures close in projection to the lens, and at 
least one of those structures (in the field of MG0751) significantly affects the lens potential. 
Subject headings: gravitational lensing — (galaxies:) quasars: individual (MG 0751+2716, BRI 
0952-0115, PG 1115+080, B1422+231, MG 1654+1346, PMN J2004-1349, 
B2114+022, HE 2149-2745) — galaxies: clusters: general — galaxies: halos 



1. INTRODUCTION 

The study of strong gravitational lens systems offers 
critical constraints on the masses, shapes, evolution, 
and substructure of galaxy dark matter halos (e.g., 
Kochanck 1991; Kecton ct al. 1998; Metcaff & MadaM 
2p01; Dalai & Kochanck 2002; Rusin ct al. 200^ 
i^jigin_^,^2£hangk .2005; .Treu fc Kooomans .2004 
Ferreras et all l2005|) . on the Hubble con stant indepen- 
dent o f the local distance ladder (e.g., 'Refsdal' '1964| 
iKochan ck & Schcchter 200^, and on the dark energy 
density (e.g.. Turner 1990; Kochanck 1996at IChaell200a: 
iLindeii 12004; .Mitchell et al., 2005 ). However, our under- 
standing of observed lenses is limited by uncertainties 
and biases in the lens models necessary to analyze the 
data. Despite improving data for lensed images and 

^This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 



lens galaxies, astrophysical applications of lensing are 
still hindered by poor knowledge of the environments in 
which strong lens systems reside. 

Several arguments suggest that lenses have complex 
environments. Statistical arguments based on galaxy 
demographics imply that at least 25% of lens galax- 
ies lie in dense environm ents such as groups and clus- 
ters (" Keeton et al.l 120001) . From spectroscopic obser- 
vations, several lenses are in fact known to lie in 
groups (MG 07 51+2716, PG 1115+080 B1422+231, 
and B1608+656 :lTonrv fc Koch anekllT999l: iKundic et alJ 
[l997ab; Tonr\l ll998r iFassnacht'' et all l2004|) . and sev- 
cral others in clusters (RX J0911+0551, Q0957+561, 
HST 14113+5221, a n d MG 2016+ 112; Knc ib ct al. 200(i 
Young et a,lJ IT9Sm iFischer et al. 1998; ISoucail et alJ 
200 1|) . Indirect evidence for the existence of other groups 
comes from the large tidal shears required to explain the 
image configurations of many four-image (quad) lenses. 



2 



Momcheva et al. 



which presumably come from mass structures near th e 
lens galaxy or along the line of sight l|Keeton et al.ll99'^ . 
The range of required shears in quad lenses could re- 
flect a range of environment densities, running from 
poor groups to rich clusters. Comparisons of the lens- 
ing rate in different surveys have also been cited as ev- 
idence that many lens galaxies probably lie in groups 
(jBlandford et al.ll2001j) . Finally, theoretical models pre- 
dict that lens galaxies reside in complex environments 
that produce substantial shears, although it is not yet 
clear whether the models pr edict shears large enough 
to explain real quad lenses ijHolder fc Schechteil 1200,'^ 
IDalal Watso;]l2nnl. 

If not handled properly, complex environments can in- 
ject uncertainties and biases into the astrophysical quan- 
tities derived from lens models (see Kccton & Zablu doS 
12004 hereafter KZ04). For example, neglecting environ- 
ment altogether leads to lens models that, for most pur- 
poses,^ are simply wrong. Approximating environmental 
effects with a simple shear term leads to models that are 
better but still tend to overestimate the Hubble constant, 
the velocity dispersion of the lens galaxy, and the dark 
energy density J^a, and to underestimate the magnifica- 
tions of the lensed images. In principle, modeling the 
full richness of environmental effects can remove these 
biases, and may also resolve the long-standing puzzle of 
why quad lenses are almost as comm on as doubles in 
statistically compl ete lens samples (seelKing et al.lll996t 
[Kpchanek 1996b; Keeton et al.lll997t iRusin fc Teefmarkl 
l200lHCohn fc Koch anek"2004'. KZ041 Such an analysis 
requires detailed knowledge of the galaxy populations, 
velocity dispersions, and projected spatial centroids of 
groups and clusters around lenses in order to determine 
how the environments affect the lens potentials. To 
date, such observations have mainly been carried out 
for the few lenses that reside in clusters, which leaves 
many lenses whose environments are known poorly or 
not at all. Worse, existing observations cannot charac- 
terize the distribution of lens environments, so we can- 
not assess environment-related biases in statistical quan- 
tities (such as 51a or the quad/double ratio) or ensemble 
properties (such as evolution or substructure). While 
the environment distribution can be predicted from the- 
oretic al models llKeetoii et al.lll997t iHolder fc Schechtieil 
120031 iDalal fc Watsonll2004|) . disagreements among the 
models, and discrepancies between the predicted distri- 
butions and the shears required to fit observed lenses, 
raise questions about the predictions. 

These issues have not been adequately addressed 
with observations, because no systematic survey of 
lens environments exists. Surveys of a few lens 
fields have been published individually or in pairs 
llYourig et alJl981bt.Kundic et al. 1997a.bl:lFischer et a~ 
199S|lTonrvll1998HTonrv Koclianeklll 999HKneib et al. 
20001 iFass nacht fc LubinI 120021 iFassnacht et alJ 120041 
Soucail et al. 2001). In many cases, though, those sur- 
veys only spanned a ~30" field around each lens, so they 
did not adequately sample group or cluster membership 
out to the virial radius (^ 0.7 Mpc for groups, corre- 

^ The important exception is measurements of the total mass 
within the Einstein radius, which ar e largely independent of as- 
surn^ tions built into lens models fe.g.. lKochanekll99UICohn et all 



sponding to ~3' at the redshifts of the lenses we study). 
We have undertaken a systematic deep and wide-field 
survey of lens fields, and here we present results for the 
first eight systems that we have targeted for multi-object 
spectroscopy. We characterize the environment within a 
~ 6' diameter field around each of the eight lenses, and 
quantify how those environments affect the lens poten- 
tials. 

Going beyond the lenses' immediate environments, we 
also consider the degree to which massive structures 
along the line of sight to a lens affect the lens potential. 
The prevalence and importance of interloping structures 
in lens fields is poorly understood. Observationally, there 
appear to be bound grou ps along the lines of sig ht to 
B0712+472 (10 members; Fassnacht & Lubin 2002) and 
MG 1131+0456 (3 members: Tonrv fc Ko chanck 2000). 
Overdensities of galaxi es are seen in the fields of sev- 
eral o ther strong lenses ()Faure et al.ll2004l iMorgan et alJ 
12005ft . but it is not yet known whether they indicate 
massive bound structures, and whether any such struc- 
tures are associated with the lens galaxies or lie elsewhere 
along the lines of sight. On the theoretical side, stud- 
ies have yielded confiicting results as to whether line- 
of-sight stru ctures are very important or negligible for 
lensing fe.g.. 'Seliak'1994': •Bar-Kana'' 1996l iKeeton et alJ 
1_997; Premadi fc Martcl 2004). Wc show here that lenses 
are, in principle, sensitive to structures over a wide range 
of redshifts and projected spatial offsets, so the practi- 
cal importance of interloping structures depends only on 
how common they are. Our photometric/spectroscopic 
pencil-beam survey of lens fields enables us to self- 
consistently identify any prominent structures at all rel- 
evant redshifts, and to assess their actual contributions 
to observed lenses. 

Separate from lensing, an important by-product of our 
survey is a sizable sample of poor groups at intermediate 
reds hifts. Only a few such samples a,re pre s ently known 
liCarlberg et alJ l2?mt IWilman et aP I2?ffl5l: IGerke et alJ 
I2005D . Groups are important laboratories for studies of 
galaxy evolution because they are the most common en- 
vironments for galaxies, and are also relatively simple 
systems in which the range of mechanisms thought to 
drive galaxy evolution (primarily galaxy-galaxy interac- 
tions) is much narrower than in hotter, denser clusters 
ijZabludoff fc Mulchaevlll998al hereafter ZM98). Unfor- 
tunately, poor groups are notoriously difficult to identify 
using conventional methods for finding clusters, due to 
their low projected surface densities, faint X-ray lumi- 
nosities, and inefficiency for weak lensing. The veloc- 
ity dispersions of nearby groups range from cr^ 200 
km for systems that are X-ray faint, late-type dom- 
inated, dynamically young, and generally similar to the 
Local Group; to ar ^ 300-500 kms^^ for systems that 
are richer. X-ray luminous, early-type dominated, and 
dynamically more evolved; up to ar ~ 1000 km s~^ for 
rich clusters. In nearby X-ray luminous groups with 
(Jr ~ 300-500 kms~^, there is always a giant elliptical 
that lies at the center of the group potential, which sug- 
gests that such galaxies form in groups via interactions 
prior to being accreted by rich clusters (ZM98). Groups 
at intermediate redshifts like those we describe here will, 
in conjunction with nearby group samples, permit us to 
observe the evolution of groups directly. 

The organization of this paper is as follows. In §2 we 
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describe our sample of eight lens systems and summarize 
previous work on them. In §3 we present our spectro- 
scopic data in the eight lens fields. In §4 we determine 
the membership, kinematics, and centroids of the groups, 
and use those properties to quantify how the environ- 
ments affect the lens models. We also explore the effects 
of line-of-sight structures on the lens models. We sum- 
marize our results and conclusions in §5. We present the 
formalism for computing the convergence and shear aris- 
ing from perturbing structures anywhere along the line 
of sight in an Appendix. Where necessary, we assume a 
cosmology with Hm — 0.3, fl\ ~ 0.7, and Hq — 70 km 



Mpc" 



2. THE SAMPLE 



Our sample consists of eight known gravitational lens 
systems with lens galaxies at intermediate redshifts 
0.25 < zi < 0.5. Four of the lenses (MG 0751-^2716, 
PG 1115-f-080, B1422-h231, and MG 1654-^1346) were 
suspected from previous studies to have complex environ- 
ments. We chose the other four lenses (BRI 0952—0115, 
PMN J2004-1349, B2114-h022, and HE 2149-2745) be- 
cause of the availability of prior imaging and photometry 
as well as accessibility from Las Campanas Observatory. 
In the remainder of this section we briefly review prior 
studies of these eight lenses. The data from this section 
are summarized in Tabled 

MG 0751-^2716 (hereafter MG0751), discovered as a 
part of the MIT-Greenbank-VLA search for gravita- 
tional lenses, is a radio lens with four images and a 
partial ring (L ehar et aL .1993,L Optical imaging of the 
system bv iLehdr et alJ (|199<j) identified an R — 21.3 
galaxy (G3) located 0.2" northeast of the brightest ra- 
dio spot as the likely lens galaxy. G3 is a satellite 
of a m uch brighter R = 1 9.1 ga laxy (Gl) located 6" 
away. iTonrv fc Kochaiiei3 l)1999tl determined the red- 
shifts of the galaxies to be zgi — 0.3501 ± 0.0003 and 
ZG3 = 0.3502 ± 0.0003. They also found a nearby emis- 
sion line galaxy to have redshift 0.3505 ± 0.0003, indi- 
cating that the lens galaxy lies in a small group with at 
least three members. Lens models by Lchar et al. (199"^ 
suggest that MG0751 requires more external shear that 
can be accounted for by the observed galaxies, which 
is consistent with the hypothesis that the lens environ- 
ment is complex. The redshift of the source quasar is 
Zs = 3.200 ± 0.001 (Tonrv & Kochanck 1999). 

BRI 0952-0115 ( hereafter BRI0952) was discovered by 
iMcMahon fc IrwinI l(T99l as a doubly imaged Zs = 4.5 
optical quasa r. The quasar is also detected at mil l imeter 
wavelengths (jOmont et al.l 11996(1 . iKeeton et all l|1998fl 
fou nd that the lens is a flattened early type galaxy, 
and lKochanek et al.l 1(20001) estimated a lens redshift of 
zi = 0.41 ± 0.05 based on fundamental plane fitting. Be- 
cause the separation between the images i s small (0.9") 
and the lens galaxy is faint (21.9 in F675W: lKeeton et al.l 
^98) , the lens redshift has not been determined spectro- 
scopically. 

PG 111 5-^080 (hereafter PG111 5) is a lens system dis- 
covered bv lWevmann et all 1(19801) , in which a radio-quiet 
quasar at redshift Zs = 1.722 is lensed into four im- 
ages (Hege^^^^^ l). The len s galaxy was first de- 
tected by Henrv & Heaslev (1986); its redshift was esti- 
mated bv iAngonin-Willaime et al.. Q993). and later im- 
proved bv lKundic et all l(1997a(l 'aiid lTonrvl l(1998(l to zi = 



0.3098±0.0002. lYoung et all l(1981a(l suggested the pres- 
ence of a smal l group of gala xies near the lens. This was 
confirmed by iKundic et all l(1997a() and iTonrvl l(1998() . 

who measured the redshifts o f a total of four gala xies 
within 20" of the lens galaxy. iK^mdic et alJ l(1997al) es- 
timated a group velocity dispersion of ar = 270 ± 70 
km s^^ from four galaxies, while Ton^ l|1998() estimated 
ar — 32 6 km s~^ from a slightly different set of four 
galaxies. iGrant et alJ 1(2004(1 detected diffuse X-ray emis- 
sion that is associated with the group and that has a 
temperature kT ~ 0.8 ± 0.2 keV; this value is consis- 
tent with typical values for low-redshift poor groups, but 
somewhat high given the measured group velocity disper- 
sion and the local (Jr-Tx relation ( Mulchaev & Zabludo^ 
119981 hereafter MZ98). PG1115 is one of nine known 
strong lens systems for which the time delay between 
different images h as been measured, so it can be used 
to determine i?o- iSchechte r et all l(1997() measured the 
light curves of the different imag es and estimated the 
time delays, and iBar-Kanal 1(1997!) then reanalyzed the 
data to give more precise results: the delay between im- 
ages B and C \s tsc = 25.0t3 g days, and the ratio of 
the delays between A (actually a combination of the close 
images Ai and A2), B, and C is IacI^ba = l-13!to.i7- 
PG1115 is one of the lenses with "anomalous" flux ratios 
thought to indicate so me sort of small-scale st r ucture 
in the lens galaxv Ce.g.. |Metcalf fc Madaull200lHGhibal 
■ c"2002': ' 



I2OOI balal fc K ochanek"2002': ' Keeton et ai]l2005a^ . 

B1422-h231 (hereafter B142 2) is a four-image lens 
discovered by IPatnaik et alJ l(1992() while searching 
for small-separation lenses among flat spectrum ra- 
dio sources i n the Jodrell Bank-VLA Astrometric Sur- 
vey fJVA S: IPatnaik et all Il992at iBrowne et all IT993 
l^ilkinson et alJ Il998t iKing et all 119991) . The source 
is a radio loud quasar at Zg = 3.62 fPa tnaik et alJ 
119921) . and the lens is a luminous ellip tical at zi = 
0.3374 (llmnev et al.lll996l IKundic etall Il997bl) . The 
lens galaxy and five nearby galaxies form a group 
at Zg — 0.338 with a rest- frame line- of-sight veloc- 
ity dis persion of = 550 ± 50 km s~^ l(Kundic et alJ 
ll997hD . Lens models for B1422 require a signifi- 
cant shear 7 ~ 0.20-0.26, w hich may be attributable 
to the group environment (Hogg fc Blandford' 119941: 
IKeeton et aIT IT997: Doblcr fc Keeton 2005). Indeed, 
from thei r estimate of t h e group 's velocity dispersion and 
centroid, iK^mdic et alJ l(1997bD estimated 7 = 0.23 and 
pointed out that the group will also create some conver- 

fice K that may affect the lens potential. IGrant et alJ 
)04) detected B1422 in X-rays (0.5-2 keV) and de- 
termined a temperature of kT — 1.0;to°3 keV, which 
is consistent with the value expected for a poor group 
(MZ98 ). B1422 is another lens with "ano malous" flux 
ratios (iMao fc Schneideilll998H Ghiba'2002D. 

MG 1654-^1346 (hereafter MG1654) was origi- 
nally detected in the M IT-Greenbank-VLA survey. 
iLangston et al.l l(1988llT989() recognized its unusual struc- 
ture in a VLA snapshot and obtained radio and optical 
mapping. The source is a = 1.74 radio quasar with a 
compact core and two extended radio lobes. The south- 
west lobe is lens ed into a ring by a 2 , } = 0.254 giant 
eUipt ic al galaxy iLLangsto n et al.lll988t iKochanek et alJ 
|2Q03). ILangston et al.l l(1989() noted an enhancement of 
the number density of galaxies near the lens; some of the 
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TABLE 1 
Gravitational Lens Galaxies 



Lens 








lb 


7 b 




Images'' ' ' 


kTx" 


iVgrp 






(J2000) 




[mag] 




[days] 




[kcV] 




[kms-1] 


MG 0751+2716 


07:51:41.46 


+27:16:31.4 


0.349=* 


21.26 


3.20 




R 




2 




BRI 0952-0115 


09:55:00.01 


-01:30:05.0 


{0.41)<: 


21.21 


4.50 




2 








PC 1115+080 


11:18:17.00 


+07:45:57.7 


0.31'' 


18.92 


1.72 


25.0+2.0 


4 


0.8+0.2 


4 


270+70 


B1422+231 


14:24:38.09 


+22:56:00.6 


0.34'^ 


19.66 


3.62 
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550+50 


MG 1654+1346 


16:54:41.83 


+13:46:22.0 


0.254=* 


17.9 


1.74 




R 






PMN J2004-1349 


20:04:07.07 


-13:49:30.7 










2 








B2114+022 


21:16:50.75 


+02:25:46.9 


0.316=* 
0.59'" 


18.63 






2+2 








HE 2149-2745 


21:52:07.44 


-27:31:50.2 


0.50'" 


19.56 


2.03 


103.0+12.0 


2 









=*Data from this work 

''Data from references in the text and from the CASTLES website jhttp : //cf a-www ■ harvard ■ edu/castlesZ 

'^Photometric redshift (Kochanek et al. 2000) 

''Number of previously known group members in addition to the lens galaxy. 
References in text. 
'R means an Einstein ring. 



nearby galaxies are comparable in brightness to the lens 
galaxy, suggesting a complex environment. 

PMN J2004-1349 (hereafter PMN2004) is a two-image 
lens disc overed in a search for radio lenses in the south- 
ern sky l)Winn et al.l I2fl0ll) . The radio spectral index 
of the images is typical for radio-loud quasars, so the 
source is considered to be a quasar despite t he lack of an 
optic al spectrum and a measu red redshi ft ("Winn et al.l 
120011) . Based on photometry, iWinn et al. '(200ijl sug- 
gested a lens redshift in the range 0.5 < zi < 1.0. Hig her- 
resolution imaging bv lWinn. Hall fc SchechteJ l)2003l) re- 
vealed a spiral lens galaxy (only the fifth one known) and 
showed that the color differences between the two im- 
ages at optical and near-infrared wavelengths can be ex- 
plained by differential extinction. The extinction analy- 
sis can be used to infer the lens redshift; it seems to imply 
somewhat low values (0.03 zi < 0.36), but that result 
depends on assumptions about the extinction curve. 

B2114+022 (hereafter B2114) was discovered a s part 
of the search for lenses in JVAS ijKing et alJll999() . Ra- 
dio maps show four distinct components within 2.4" in 
a configuration that is atypical for lenses. Furthermore, 
the sources can be divided into two pairs with distinct 
radio surface brightnesses and radio spectra: sources A 
and D are similar to each other, and sources B and C 
are similar to one another, but the two pairs are clearly 
different. Ground-based and HST optical imaging and 
spectroscopy do not detect the lensed images but re- 
veal two le ns galaxies at zn — 0.3157 and zi2 = 0.5883 
ijAugustcTe t al. 20 0ll), suggesting a c omplex lensing ge- 
ometrv. IChae. JMaofc Augustd 1)20011) could explain two 
of the radio components (A and D) as lensed images us- 
ing a two-plane lens model. It is not known whether 
the other components (B and C) are images of the same 
source (unlikely), images of a different source, or struc- 
ture related to the Gl lens galaxy. No lens models that 
explain these components have been published. 

HE 2149-2745 (hereafter HE2149) is a doubly im- 
aged broad absorption line qu asar at redshift z,^ = ^ 
2.033, which was discovered by iWisotzki et alJ ij 19961) 
in the Hamburg/ESO wide-angle survey for bright 



quasars. iBurud et alJ l)2002t) reported a redshift of 
zi = 0.495 ± 0.01 for the elliptical lens galaxy, consistent 
with the photometric r e dshift estimat e of zi = 0.43tQ;gg 
ijKochanek et al.ll2000|) . IBurud et alJ also measured the 
time delay between the two images to be At — 103 ± 12 
days. Based on the large number of red non-stellar 
objects in R-band images of the field around the lens, 
iLopez. Wucknitz fc Wisotzkil fl99^ suggested that the 
lens galaxy might be a member of a cluster. 

3. SPECTROSCOPIC DATA 

We first identified galaxies for follow-up spectroscopy 
from two-color, wide-field imaging of each lens field. We 
obtained deep images in I and either V or R during the 
period from IVIay 2002 to June 2004 using the 36' x 36' 
Mosaic Imager on the 4-m telescopes at Kitt Peak Na- 
tional Observatory and Cerro Tololo Inter- American Ob- 
servatory. We reduced these images and extracted photo- 
metric catalogs followin g standard methods using IRAF^ 
and SExtractor (Berlin fc Arnouts 1996). A more de- 
tailed description of the photometric analysis is presented 
in a separate paper ijWiniams et alj[2005i1 . 

We selected spectroscopy targets using a prioritization 
scheme based on objects' colors and projected distances 
from the lens. Highest priority was given to objects 
populating a red sequence in the color-magnitude dia- 
gram that is consistent with the lens redshift, and to 
targets that lie within a group-like virial radius of 0.7 
Mpc (ZM98), which corresponds to ^ 3' over the red- 
shift range of our lenses. We applied a magnitude cut 
at / = 21.5 to assure reasonable exposure times. This 
limiting magnitude corresponds to I* + 4.2 at the low 
redshift limit of our sample (z = 0.25) and to /* + 2.8 at 
the high redshift limit (z = 0.5), where /* is the observed 
magnitude of an galaxy, adopted from Williams et ^ 
(^005). We obtained multislit spectroscopy during two 
observing runs, March 1-4 and August 30-September 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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2, 2003, with the Low Dispersion Su rvey Spectrograph 
rLDSS-2: fATrington-Smith et alJlflflflD at the 6.5-m Mag- 
ehan 2 (Clay) telescope at Las Campanas Observatory. 
All spectra were taken with the medium blue grism (300 
1/mm; 5000 A blaze) over a wavelength range of 3900- 
8OOOA. We used 1.03" slitlets, resulting in a spectral res- 
olution of I5A FWHM. Each slitmask had dimensions 
~ 5' X 7', typically contained 20-30 targets, and was ob- 
served for 4x900 s. 

Our sky coverage is shown in Figure ^ Each panel 
spans 15' x 15' and is centered on the lens galaxy. Rect- 
angles show the outlines of our slitmasks. There is a 
noticeable difference in the sampling between the two ob- 
serving runs. In March 2003, when we observed MG0751, 
BRI0952, PG1115, B1422, and MG1654, the masks lay 
mostly on top of one another, providing exhaustive cov- 
erage of the immediate surroundings of each lens galaxy. 
Because most of our fields were far north for Magellan, 
we aligned the masks with the paralactic angle to min- 
imize the effects of atmospheric dispersion. In contrast, 
the masks for MG1654, PMN2004, B2114, and HE2149 
in the August 2003 observing run were tiled to cover a 
larger area of the sky but still overlap significantly in 
the 3' projected radius around the lens. No attempt was 
made to align the masks with the paralactic angle be- 
cause the targets were close to or south of the celestial 
equator and were therefore observed at relatively low air- 
mass. 

The figure shows that our sampling is very good within 
3' around the lens for all fields, with a minor excep- 
tion in B1422. We observed four masks for each of 
MG0751, BRI0952, PG1115 and B1422, five masks for 
each of PMN2004, B2114 and HE2149, and seven masks 
for MG1654 (which was observed during both runs). We 
discuss our spectroscopic completeness below. 

We reduced all spectra using standard IRAF proce- 
dures and corrected them to the local standard of rest 
using the IRAF routines RVCORRECT and DOPCOR. 
We determined the radial velocities using the cross- 
correlation of absorption lines (XCSAO) and/or using 
emission line identifications (EMSAO)^. If both emis- 
sion and absorption line velocities were found, the quoted 
value is a weighted average of the two. Marc Postman 
kindl y provided template galaxy spectra ( Postma n et al.l 
I2OO2I and private communications). We visually in- 
spected every fit to ensure accuracy. 

The number of objects observed in each field is listed 
in Table 121 Galaxies represent ^42% of the objects tar- 
geted spectroscopically, while stars originally misclassi- 
fied as galaxies are '^22%. "Failed" targets (~32%) are 
those for which were unable to obtain velocities. The 
majority of failed targets, especially at fainter magni- 
tudes, were absorption-line systems for which the signal- 
to-noise was too low to allow successful cross-correlation. 
Other causes of failures are low surface brightness or poor 
astrometry (more problematic in March than August). 
The fraction of stars was significantly lower in August 
(15%) than in March (26%), thanks to improvements 
in s tar-galaxy separation made between the two runs 
fsee 'Williams et 311120051 ). The large number of stars 
in PMN2004 is due to the high stellar density at this 

^ XCSAO and EMS AO are routines in the RVSAO IRAF package 
IKurtz fc MinldlTOQl) 



relatively low galactic latitude and longitude {I = 28°, 
b= -22°). 

Figure |21 shows the completeness of our spectroscopy 
with respect to the photometry. The solid-line histogram 
shows the magnitude distribution for all galaxies in the 
photometric catalog projected within 3' of the lens, while 
the shaded histogram shows the subset for which we de- 
termined velocities. Our target selection scheme, based 
on colors and projected offsets from the lens, misses a 
few of the brightest galaxies (many of which are likely to 
be foreground objects). We miss a larger fraction of ob- 
jects at the faint end, but these are very few and/or less 
massive galaxies, which would significantly affect the lens 
potential only if very close to the lens galaxy. In total, 
there are only four galaxies within 10" of the lens galax- 
ies (the zone in which a small pcrturber can have even a 
moderate effect; see Fig.lHlbelow) that are present in our 
photometric catalog but have no determined velocities. 
Two of those (in B2114 and PMN2004) are below our 
spectroscopic magnitude limit of / = 21.5. The other 
two galaxies are in the field of IIE2149, have I magni- 
tudes of 19.97 and 21.08, and lie 7" and 10" respectively 
from the lens. In color-magnitude space l) Williams et alJ 
12005(1 they lie on red sequences identified as line-of-sight 
structures at z = 0.45 and z = 0.60 (see §4.5). 

We were never able to put slits on all of our highest- 
priority objects; this limitation is inherent to multislit 
spectroscopy. For lensing purposes, the main effect of 
spectroscopic incompleteness is to cause us to underesti- 
mate environment-related lensing biases (see §§4.4-4.5). 
In particular, the fact that we prioritize galaxies thought 
to lie at the lens redshift (color selection) means that we 
undersample line-of-sight structures and hence underes- 
timate their contributions to the lens potential. 

As already mentioned, we improved our photometric 
catalogs after some of the spectroscopic targets were se- 
lected. As a result, there are 15 galaxies whose velocities 
we measured that do not actually appear in the final pho- 
tometric catalog. (These galaxies lie under bleed trails in 
the imaging, so accurate photometry is not possible.) In 
addition, there are another 25 galaxies whose velocities 
we measured that are (mis-) classified in the final pho- 
tometric catalog as stars. We omit all of these galaxies 
when comparing the spectroscopic and photometric cat- 
alogs for the purpose of understanding our spectroscopic 
completeness (i.e., they are excluded from Fig. How- 
ever, these remain a part of our spectroscopic catalog 
and all analyses based on that catalog. 

To estimate the zero-point velocity correction and ex- 
ternal velocity errors, we cross-correlate 403 sky spectra 
extracted from our data with the same templates used for 
the galaxy spectra and find a mean velocity of u = 40±50 
km s~^. We also determine the velocities of 153 of the 
serendipitously observed stars and find v = 30 ± 180 
km s~^. Both methods give mean velocities comparable 
to or smaller than the dispersion, and much smaller than 
the velocities of the objects in the sample. We therefore 
conclude that no zero point correction is needed. 

Table 121 lists our spectroscopic catalog. For each entry 
we give the catalog name, J2000 coordinates calibrated 
to USNO-B2.0, projected distance from the lens in ar- 
cmin, aperture magnitude within a fixed physical size of 
~6.5 kpc fsee Wil liams et al 2005t') . the heliocentric ra- 
dial velocity, velocity error, redshift and redshift error. 
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Fig. 1. — Sky coverage for LDSS-2 multislit spectroscopy in the fields of the eight lens systems in our sample. Each panel is 15' X 15', 
centered on the lens galaxy. Rectangles show the positions of our slit masks. Each slitmask covers ~ 5' X 7' and includes 20-30 slitlets 
parallel to the long side of the mask. The 3' circle around the lens galaxy corresponds to ~ 0.7 Mpc, a group-like virial at the redshifs of 
the lenses in our sample. We observed four masks for each of MG0751, BRI0952, PG1115, and B1422, five masks for each of PMN2004, 
B2114, and HE2149, and seven masks for MG1654. The masks from the first observing run lack a particular sampling pattern but provide 
excellent coverage of the expected group virial radius in all cases except B1422, where a small portion remains unsampled. The masks 
from the second observing run were tiled in an attempt to maximize the sky coverage while still providing good sampling of the immediate 
surroundings of the lens. 



TABLE 2 
LDSS-2 Observations 



Lens 


Date 


# Galaxies^ 


# Stars 


# Failed 


Total 


MG0751 


March 2003 


38 


24 


27 


89 


BRI0952 


March 2003 


47 


24 


19 


90 


PG1115 


March 2003 


47 


28 


11 


86 


B1422 


March 2003 


53 


14 


26 


93 


MG1654 


March 2003 


39 


30 


32 


101 




August 2003 


20 


6 


35 


61 


PMN2004 


August 2003 


41 


43 


35 


119 


B2114 


August 2003 


38 


1 


46 


85 


HE2149 


August 2003 


41 


8 


66 


115 



^Including QSOs. 



The last column describes the method from which the 
velocity was obtained: 1 for absorption lines, 2 for emis- 
sion lines, or 3 for a combination of both. Missing data 
means that the object was not in the final photometric 
catalog as described above. Such objects have no iden- 
tification names and are numbered successively starting 
with 90001. 

4. RESULTS AND DISCUSSION 

4.1. Group Membership 

The environments of most strong lenses are not well 
characterized. Our first goal is to determine whether 
each lens galaxy lies in a group or cluster, and if so to 
identify the other member galaxies. Even in cases where 
groups were already identified (MG0751, PG1115, and 



B1422), the number of group members known previously 
ranged from three to six. Our deep, wide-field spectro- 
scopic sampling has the potential not only to find new 
groups, but also to increase the membership of known 
lens groups to the point where robust determinations of 
the group velocity dispersions and centroids are possi- 
ble. These are essential for understanding how a group 
affects the lens potential, as discussed in §4.4. Addition- 
ally, a more complete inventory of the brightest (most 
massive) group members, and their contributions to the 
lens potential, will also greatly improve lens models. 

We present the redshift histograms for all eight of our 
lens fields in Figure 13 In the left panels, the shaded his- 
tograms include all galaxies that lie within a projected ra- 
dius of 3' (a group-like virial radius) around the lens, and 
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Fig. 2. — Apparent magnitude histograms of galaxies in the photometric catalog (solid line), and of those galaxies within 3' of the lens 
for which redshifts were obtained (shaded). Our target selection scheme, based on colors and projected offsets from the lens, misses a few 
of the brightest galaxies. We miss a larger fraction of objects at the faint end, but these are presumably less massive galaxies that do not 
contribute significantly to the lens potential. 



TABLE 3 

Lens Field Galaxy Properties — EXAMPLE 



ID 


a 

[hh:mm;ss] 


5 

[dd:mm:ss 




b 
['] 




I 

[mag 




cz 

[km s~ 


Acz 


z 




Spectral 
Type 


MG0751: MWKZ 


GAL 






















... 9809 


07:51:30.20 


27:14:43 


1 


3 


10 


20 


9 


168070 


30 


0.56023 


0.000083 


2 


... 9739 


07:51:30.77 


27:17:55 


7 


2 


76 


20 


6 


105970 


70 


0.35323 


0.000237 


2 


... 9726 


07:51:31.00 


27:14:16 


7 


3 


24 


20 


5 


175630 


50 


0.58543 


0.000173 


3 


... 9582 


07:51:32.04 


27:12:59 


2 


3 


30 


19 





91080 


30 


0.30360 


0.000113 


3 


... 9606 


07:51:32.28 


27:12:58 


9 


4 


10 


20 


8 


106950 


100 


0.35650 


0.000320 


1 


... 9570 


07:51:32.49 


27:17:38 


6 


2 


29 


20 


2 


104680 


30 


0.34893 


0.000113 


3 


... 9560 


07:51:32.73 


27:14:42 


1 


2 


67 


21 





106300 


30 


0.35433 


0.000113 


2 


... 9333 


07:51:34.61 


27:15:45 


5 


1 


71 


19 


5 


91220 


50 


0.30407 


0.000153 


1 


... 9297 


07:51:34.66 


27:13:37 


7 


3 


28 


18 


4 


28910 


40 


0.09637 


0.000133 


1 


... 9274 


07:51:35.24 


27:17:37 


1 


1 


76 


19 


6 


105350 


110 


0.35117 


0.000380 


1 


... 9225 


07:51:35.43 


27:17:07 


9 


1 


48 


19 





79820 


60 


0.26607 


0.000213 


3 


... 9239 


07:51:35.75 


27:15:22 


8 


1 


72 


20 


9 


104960 


80 


0.34987 


0.000280 


1 


... 9100 


07:51:36.66 


27:19:39 


8 


3 


31 


18 


4 


74510 


50 


0.24837 


0.000160 


1 


... 9047 


07:51:36.96 


27:19:27 


9 


3 


10 


18 


3 


104220 


40 


0.34740 


0.000140 


2 


... 9120 


07:51:36.98 


27:18:40 


6 


2 


37 


20 


9 


72100 


60 


0.24033 


0.000203 


2 


... 9049 


07:51:37.73 


27:18:06 


5 


1 


78 


20 


4 


71940 


50 


0.23980 


0.000163 


2 


... 8794 


07:51:37.92 


27:16:12 


9 





85 








167590 


100 


0.55863 


0.000320 


1 


... 9006 


07:51:38.04 


27:17:33 


7 


1 


28 


20 


1 


168620 


90 


0.56207 


0.000307 


1 


... 8816 


07:51:40.32 


27:16:22 


1 





31 


20 


9 


104940 


100 


0.34980 


0.000323 


1 


... 8669 


07:51:41.07 


27:19:43 


7 


3 


20 


18 


9 


147320 


120 


0.49107 


0.000390 


1 


... 8682 


07:51:41.50 


27:16:31 


9 





00 


20 


1 


104810 


120 


0.34937 


0.000410 


1 


... 8673 


07:51:41.84 


27:16:29 


2 





09 


21 


5 


167090 


10 


0.55697 


0.000037 


2 


... 8476 


07:51:43.25 


27:17:53 


3 


1 


41 


19 


8 


169470 


210 


0.56490 


0.000687 


1 


... 7921 


07:51:43.26 


27:16:06 


3 





58 








124780 


30 


0.41593 


0.000103 


3 


... 8385 


07:51:44.18 


27:16:39 


7 





61 


20 


7 


112430 


200 


0.37477 


0.000657 


1 


... 8257 


07:51:45.32 


27:17:02 


4 





99 


19 


4 


105110 


60 


0.35037 


0.000200 


1 


... 8288 


07:51:45.50 


27:18:51 


4 


2 


49 


20 


4 


60380 


100 


0.20127 


0.000317 


1 



Note. — This table is published in its entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. Objects with missing data were not found in our final 
photometric catalog for reasons explained in the text. 
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A'tot indicates the number of galaxies in the histogram. 
This number excludes a few high-redshift AGNs that fall 
outside the range of the plot. The histograms include 
all galaxies with spectroscopic redshifts both from our 
sample and from the literature. In particular, we have 
added the following 16 galaxies to our catalog : Gl, G6, 
and G7 in MG0751 ijTonrv fc Kochanekl 11991: the lens 
galaxy GL as well as Gx in PG1115 l)Tonrvlll998|) : the 
lens galaxy Gl as well as G2. G3, G4, G6, G8, G9, GIO, 
and Gx in B1422 l|Kundic et al| 1997b: Tonrv"1998'): the 
zi = 0.59 lens galaxy in B2114 (|Augusto et al. 2001); and 
the z i = 0.495±0.01 lens galaxy in HE2149 (Burud et al] 
I20n2|) . We have not included the lens galaxy in BRI0952 
(because it only has a photometric redshift estimate) or 
in PMN2004 (because no good redshift estimate exists). 
The vertical line shows the position of the lens galaxy 
as listed in Tabled In six cases (see Table QJ, there is 
clearly a peak in redshift space at or near the lens galaxy 
redshift. 

We determine the group membership by applying 
a pessimist ic Scr clipping algorithm (as suggested by 
lYahil fc Vid al 197'^ to any redshift peak containing a 
lens galaxy. This procedure removes the galaxy most de- 
viant from the mean redshift and recalculates the mean 
and velocity dispersion of the distribution. If the omit- 
ted galaxy is more than three new standard deviations 
away from the recomputed mean, it is rejected. This loop 
is executed until an omitted galaxy is not rejected. We 
use statistical bi-weight estimators, which are more ro- 
bust for small sample sizes than the standard estimators 
(p3ccrs ct al. 1990), to calculate the location (mean red- 
shift) and scale (velocity dispersion). In Figure A^g^p 
is the number of group members determined by the 3cr 
clipping algorithm. 

The right-hand panels in Figure 13 show a cut within 
±5000 km s^^ of the lens galaxy velocity to present a bet- 
ter view of the lens groups themselves. The peak height 
is lower than in the full histogram because the bin size 
is smaller, but groups are still easily recognizable in six 
cases: MG0751, BRI0952, PG1115, B1422, MG1654 and 
B2114. The groups m BRI0952, MG1654, and B2114 
are new discoveries. For the previously known groups, 
we have increased the number of group members from 
three to 13 m MG0751, from five to 13 m PG1115, and 
from six to 16 in B1422. The projected spatial distribu- 
tions of galaxies in the six groups are shown in Figure 

H 

In the case of BRI0952, we suggest that the lens 
galaxy may belong to the five-member group at Zg = 
0.422, which is consistent with its photometric redshift 
estimate of zi = 0.41 ± 0.05. In B1422, the group 
seems to consist of two clumps — one around the lens 
gal axy and the other to the nort heast of it (see Figure 
El). iZabludoff fc MulchaevI l|1998bD have found analogous 
substructure in nearby groups. B1422 is our best sam- 
pled group (16 members) and we expect that increasing 
the membership of the other groups might reveal similar 
dumpiness. In B2114, the group is associated with the 
foreground lens galaxy (the first vertical dashed Hue). 
There is a second peak slightly in front of the group 
around the lens galaxy, but the galaxies are significantly 
offset from the lens on the sky and thus not important 
for lensing. We classify as "group members" only the 
galaxies in the peak at the lens redshift. 



Our spectroscopic findings agre e well with the ex pec- 
tations set by our photometry fsee lWilliams et al.12 005'1. 
All six lenses where we find groups show a compelling 
red sequence at the lens redshift. In IIE2149, the color- 
magnitude diagram shows a well defined red sequence 
corresponding to 2 ~ 0.28, where we see a prominent 
hne-of-sight structure in FigureEl In PMN2004, the lack 
of compelling structures in either the color-magnitude 
diagram or the redshift histogram suggests that there 
are no significant structures along the line of sight, al- 
though it is difficult to draw firm conclusions from non- 
detections (especially given incomplete spectroscopy) . In 
any case, the agreement between our photometric and 
spectroscopic results reassures us that we understand the 
data and their implications. These results are discussed 
further in iWiUiams e t al. (2005). 

4.2. Group Kinematics 

The group velocity dispersion provides a key observ- 
able probe of the group potential and its effect on 
lens models. Group velocity dispersions based on small 
member catalogs are uncertain and may be biased, be- 
cause poor sampling of the underlying velocity distribu- 
tion tends to underestimate the true velocity dispersion 
(ZM98). Our more extensive member catalogs now al- 
low us to measure the group velocity dispersions more 
accurately than was previously possible. 

To determine the mean velocities, €», and line-of-sight 
velocity dispersions, (7^, of the six groups in our sam- 
ple, we use bi-w eight estimators of location and scale 
ijBeers et al.l 11990) because of their superiority at de- 
weighting tails in the velocity distribution. In BRI0952, 
which has only five members, the bi-weight estimator 
routine fails, so we use the standard method for cal- 
culating the mean velocity and velocity dispersion from 
Danese et al. (1980). In general, standard and bi-weight 
methods yield similar means and velocity dispersions. 
We use all known group members, including the 10 found 
in the literature, to determine these kinematic properties. 
We apply a standard 1/(1 -I- Zg) cosmological correction 
to the velocity dispersions (Beers ct al. 1990). The kine- 
matic properties of the groups are presented in Table 0] 

The six groups have velocity dispersions ranging from 
not™ to 470l:go° km s-^ In the nearby universe, this 
range of velocity dispersions describes systems running 
from dynamically young, unrelaxed systems like the Lo- 
cal Group, up to more dynamically relaxed, X-ray lumi- 
nous groups. We are probably seeing a similar range of 
groups in our intermediate-redshift sample. It is impor- 
tant to note that the total number of galaxies for which 
velocities were obtained is similar in all six fields, and 
that A^grp is roughly correlated with (7^, suggesting that 
the differences between the measured ar values are real 
and not due to variable sampling. 

As a further check on the accuracy of our group ve- 
locity dispersions, we can determine the X-ray tempera- 
ture that would be derived by combining our measured 
fJr values with the Ur-Tx relation for nearby groups and 
clusters (ZM98), and compare that with the temperature 
measured directly from X-ray observations. Only two of 
the groups in our sample have been observed in X-ray s 
with Ghandra: PG1115 and B1422 ijGrant et al J 120041) . 
The X-ray temperatures expected from the ar-Tx rela- 
tion (1.5 keV and 1.7 keV, respectively) are consistent 
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TABLE 4 

Group Kinematic Properties 



Lens 




Ngrp 




C a,ce,n 

["] 




["] 


[ km 




V 

[ km s 


Sv 


z 


[ km 


5(7 r.grp 


MG0751 


39 


13 


07:51:40.7 


±11 


+27:16:53 


±9 


104220 


105970 


104980 


±100 


0.3499 


320 


+ 170 
-110 


BRI0952 


44 


5 


09:54:56.1 


±34 


-01:29:58 


±27 


125000 


127000 


126510 


±30 


0.4217 


170^ 


+ 150 
-100 


PG1115 


48 


13 


11:18:16.8 


±11 


+07:45:36 


±9 


92120 


93960 


92970 


±110 


0.3090 


440 


+90 
-80 


B1422 


57 


16 


14:24:41.0 


±11 


+22:55:42 


±9 


100640 


102810 


101540 


±130 


0.3385 


470 


+ 100 
-90 


MG1654 


59 


7 


16:54:39.3 


±27 


+13:47:15 


±20 


75390 


76210 


75750 


±100 


0.2525 


200 


+120 
-80 


B2114 


38 


5 


21:16:51.4 


±34 


+02:10:59 


±27 


93000 


95000 


94240 


±80 


0.3141 


110 


+ 170 

-80 



'^Velocity dispers ion calculated in the manner of IDanese et al.l 11980) . instead of using the statistical bi- weight estimator of scale 
IBeers et aLHTggoh . See text. 
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Fig. 3. — (Left) Galaxy redshift distributions of the eight fields in our sample. The bin size is lOOOkms-l. The shaded histo gram 
includes all galaxies that lie within a projected radius of 3' about the lens. The vertical dashed lines show spectroscopic lens galaxy 
redshifts from the literature. The vertical dotted line for BRI0952 shows a photometric estimate of the lens galaxy redshift. In PMN2004, 
the vertical dot-dashed lines show two different model-implied estimates of the lens galaxy redshift I Winn. Hall fc SchechteE.t2003 ). A^tot is 
the total number of galaxies included in the histogram, while A^grp is the total number of group members. (Right) A close-up of the range 
±5000 km s~^ centered on the mean group velocity. The bin size is 500 km . The shaded histogram shows confirmed group members. 



within the 95% confidence hmit with the observed val- 
ues (0.8 keV and 1.0 keV, respectively). This suggests 
that our values of cr^ for these groups are reasonable. 

4.3. Group Centroids 



The projected offset of the lens galaxy from the group 
centroid is another key ingredient in estimating the con- 
tribution of a group to the lens potential. The position 
of the brightest group galaxy relative to the spatial and 
kinematic centroid of the group is also an important con- 
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Fig. 131 — continued. 



straint on models of giant elliptical formation (ZM98). 
(The lens galaxy may or may not be the brightest group 
galaxy, as discussed below.) In this section we calcu- 
late the projected spatial and kinematic centroids of the 
groups by averaging the sky positions and velocities, re- 
spectively, of all group members. For members whose 
velocities were added from the literature, we use coordi- 
nates from our own photometric catalog when possible in 
order to maintain a consistent coordinate system. Nei- 
ther the projected spatial centroid nor the mean velocity 
is weighted by the luminosity, because we do not want 
to introduce an a priori bias toward the brightest group 
galaxy or assume implicitly that the mass-to-light ratios 
for all group members are the same. Nevertheless, the 
luminosity-weighted centroids are within 2(t of the un- 
weighted centroids for all groups except B1422 (where 
four bright galaxies close to the lens pull the luminosity- 
weighted centroid 3.5cr away from the unweighted cen- 
troid). 

4.3.1. Projected Offset Between the Spatial Centroid 
and the Lens Galaxy 

If a group around a lens has a common dark matter 
halo, the effects of that halo on the lens potential are 
sensitive to any projected spatial offset between the halo 



centroid and the lens galaxy. In particular, the offset de- 
termines the relative importance of convergence (gravi- 
tational focusing) and shear (tidal distortions) from the 
group halo. (See §4.4 and the Appendix for details.) 
While we obviously want to determine the offset in each 
lens/group system, we also seek to understand the distri- 
bution of offsets because that affects the distribution of 
convergence and shear, which in turn affects statistical 
applications like constraining the dark energy or under- 
standing the quad/double ratio. For MG0751, PG1115, 
and B1422, the new members we have found allow us to 
measure the offsets more precisely. For the newly discov- 
ered groups around BRI0952, MG1654, and B2114, we 
are able determine the offsets for the first time. 

Figure^ shows that in some cases there is a clear offset 
between the lens galaxy and the projected group spatial 
centroid. In B1422, the offset is substantial; the lens 
galaxy lies outside the 2a errors for the group centroid 
and is not the galaxy closest to the group centroid. In 
MG0751, PG1115 and MG1654 the lens galaxy is only 
marginally within the 2a centroid errorbars, and is also 
not the galaxy closest to the centroid. We use these 
spatial offsets in our calculations of the shear and the 
convergence due to the lens environment in §4.4. 
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Fig. 4. — Spatial distribution of the group member galaxies on the sky. North is up and east is to the left. The fields are centered on the 
lens. Each panel has an angular size of 6.8' X 6.8', which is roughly equivalent to the typical size of nearby poor groups. The open circles 
mark the group galaxies, an open square denotes the lens galaxy, and a four-pointed star indicates the group centroid and its 2cr errorbars 
(see §4.3). All galaxies with measured velocities are marked with small solid points. The scale bar in the lower left corner of each panel 
corresponds to 200 kpc at the lens redshift. 



4.3.2. Brightest Group Galaxy vs. Group Potential 

Key issues that bear on the evolution of groups and 
their galaxies are whether the brightest group galaxy 
(hereafter BGG) is kinematically and spatially distinct 
from (1) the other group members or (2) from the cen- 
ter of the group potential. In nearby X-ray luminous 
groups ((Tj. > 300 km s~^), there is a bright, giant ellip- 
tical galaxy that occupies the center of the potential (as 
defined by the spatial and kinematic centroids) and that 
lies on an orbit distinct from the other group members 
(ZM98). Our survey now allows us to ask the same ques- 
tions for intermediate-redshift groups, and to consider 
what the answers imply about group evolution. 

In each of our six groups the BGG appears to have an 
early-type morphology. The identifications are as follows: 

• MG0751: The BGG is the Gl galaxy, which hes 
6" fr om the lens galaxy and is 2.2 mag brighter 
in I. iLehar et al.l l|l997fl fit elliptical profiles to 
Gl and find an a cceptable fit. Furthermore, 
iTonrv fc Kochanekl (|T999) find that Gl is much 
brighter that the lens galaxy and has a pure ab- 
sorption line spectrum. 

• BRI0952: We identify the BGG by visual inspec- 
tion of our photometry as an elliptical 2.3' away 
from the lens galaxy and 1.5 mag brighter in I. 

• PG1115: The BGG is the giant early-type galaxy 
labeled Gl, located 12" away from the lens galaxy 



and 0.7 mag brighter in I. IImpev et alJ ljl998j) com- 
ment on its early type morphology. 

• B1422: iKundic et all (p97b) claim the BGG is the 

galaxy G3 located 8" from the lens galaxy and 1.5 
mag brighter in V. It is the brightest group galaxy 
in our sample, has a central location, and exhibits 
an early type spectrum. 

• MG1654 and B2114: In both cases, the lens galax- 
ies are the brightest galaxies in the groups and are 
also ellipticals as clas^ sified bv lKochanekl ljl995j) and 
lAueusto et al.l l)2001[) . respectively. 

The top panel of Figure |31 shows the projected spatial 
X and kinematic y offsets from the group centroid for the 
BGGs (filled squares) and for all other member galax- 
ies (filled circles) in our sample groups. The y errorbars 
(ey) represent the la uncertainties based on adding the 
galaxy velocity errors and the group mean velocity er- 
rors in quadrature. The x errorbars (ex) represent the 
la uncertainties based on the adding the centroid errors 
and the individual galaxy position errors in quadrature 
(the centroid errors dominate). To estimate the statisti- 
cal errors of the centroid for a group of A^g^p members, 
we carry out a statistical bootstrap analysis where we 
draw 500 random samples of Agrp galaxies without re- 
placement from the B1422 group (the one with the most 
members), and adopt the variance of the centroid posi- 
tion as its error. For B1422, we use the smallest of the 
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errors calculated for the other groups. 

We define the quantity — (x/Sx)^ + (y/Sy)^ as a 
measure of the phase-space distance of a galaxy from the 
group centroid. Here Sx and dy are the rms deviations in 
X and y for all galaxies plotted in the top panel of Figure 
A galaxy will have a large value of R if it has a large 
peculiar velocity and/or a position that is far from the 
projected spatial centroid. Conversely, galaxies at rest 
in the center of the group potential will have small R 
values. The bottom panel of Figure O shows the distri- 
butions of i? values for the BGGs (shaded histogram) and 
for all other group members (unshaded histogram). We 
can now use these distributions to answer two questions 
about the BGGs. 

Are the BGGs distributed differently than the other 
group galaxies? We compute the R distributions for 
BGGs and for all other group galaxies (see the bottom 
panel of Figure [SJ , and then compare them using three 
statistical tests: the KS-test (to compare the overall dis- 
tribution), the t-test (to compare the means), and the 
F-test (to compare the variances). All three fail to dis- 
tinguish between the two distributions. In other words, 
there are no significant differences between the orbits of 
BGGs and the orbits of other group members, at least 
for these small samples. This result differs from observa- 
tions of nearby groups (see ZM98). It is not clear from 
our present sample whether our result indicates real evo- 
lution in the group galaxy population between z ~ 0.3 
and z = 0, or is due simply to small number statistics. 
It will be interesting to return to this question with the 
larger sample of lensing-selected groups that we are ob- 
taining. 

Are the BGGs consistent with the group centroids? We 
compare the distribution of R values for the BGGs with 
a model distribution expected for a galaxy lying at the 
bottom of the group potential. To incorporate measure- 
ment errors, we treat the model distribution as a Gaus- 
sian in X and y with rms deviations of ex/Sx and £y/Sy, 
respectively, and make 1000 random draws using the ap- 
propriate values of ex and for each group. The bottom 
panel of FigureElshows the observed and model distribu- 
tions. A KS-test gives 6.6 x 10""^ as the probability that 
the two samples are drawn from the same distribution. 
The probability that the means of the two distributions 
are the same is 7 x 10~^, while the probability that the 
variances are the same is 10~^. We conclude that the 
BGGs do not occupy the center of the group potential. 

This offset of the BGG from the kinematic and spatial 
centroid of the group is not seen in nearby X-ray lumi- 
nous groups. While our result could suggest group evo- 
lution from z ^ 0.3 to now, another possibility is that we 
are not comparing apples to apples. For example, there 
are groups in our sample with velocity dispersions lower 
than what is typical for dynamically-evolved X-ray lu- 
minous groups nearby. Among nearby groups, lower-cr^ 
systems tend to be dynamically younger and are more 
likely to have an offset BGG (e.g., the Local Group). 
This may be true among our z 0.3 groups as well. 

Support for this latter interpretation comes from a 
closer look at the BGGs in the three high-Cr groups in 
our sample: B1422, PG1115, and MG0751. The BGGs 
of MG0751 and PG1115 are the galaxies closest to the 
projected spatial/kinematic centroid in their respective 
groups, and, within the errors, are consistent with be- 
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Fig. 5. — (Top) Projected spatial and kinematic offsets of the 
brightest group galaxy (BGG; filled squares) and all other group 
members (filled circles) from the group centroid for five of our 
groups. The lens galaxies are marked with open circles. (The 
lens galaxy in BRI0952 is omitted because it does not have a mea- 
sured redshift.) The velocity offset is normalized by the velocity 
dispersion of the group to compensate for the differences among 
the group potentials. The y errorbars represent the 68% confi- 
dence level based on adding the errors for the lens galaxy velocity 
and the mean group velocity in quadrature. The x error bars are 
the 68% confidence level based on a statistical bootstrap test (see 
text). (Bottom) Distribution of the phase-space offset R for the 
BGGs (shaded histogram) and all other group galaxies (open his- 
togram, normalized by the number of the brightest group galaxies). 
Statistical tests fail to distinguish between these two distributions. 
The heavy line (scaled down by a factor of 400 to fit the y-axis) 
shows the model R distribution for a galaxy assumed to lie at the 
bottom of the group potential (see text). A KS-test gives 7 X 10""^ 
as the probability that the BGG and model distributions are drawn 
from the same parent distribution, and a t-test gives 7 X 10^''* as 
the probability that the means of the BGG and model distributions 
are the same. We conclude that the BGGs generally do not occupy 
the center of the potential, and, that, within the large uncertain- 
ties, their phase space distribution is consistent with that of the 
other group members. 



ing at the centroid. In contrast, the BGG in B1422 lies 
at a large projected distance from the projected spatial 
centroid (^--^0.22 Mpc) and has a substantial peculiar ve- 
locity (1270 km s""'^). The BGG spatial and kinematic 
offset, together with the dumpiness of B1422 (see §4.1), 
suggest that this system is not yet relaxed. Therefore, in 
two out of three cases, the high-cTr groups in our sample 
are comparable to those nearby groups with high velocity 
dispersions and centrally located BGGs. 

We also compare the projected spatial centroid with 
the peak of emission from the diffuse, luminous X-ray 
halos in PG1115 and B1422 (the two highest-Cr groups) 
seen bv iGrant et all l|2004f) . In PG1115, the projected 
spatial centroid and X-ray peak are consistent within 2a. 
By contrast, in B1422 the peak of the X-ray emission is 
substantially offset (more than 3cr) from both the un- 
weighted and luminosity-weighted group centroids. 

Our ability to address the questions of evolution raised 
in this section will improve when we finish obtaining 
the larger sample of lensing-selected groups. Since lens- 
ing is sensitive to groups spanning the redshift range 
0-2 ^ z < 1, it will naturally provide the large red- 
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shift baseline needed to probe evolution. Working with a 
self-consistently selected sample of groups will mitigate 
selection effects. 

4.4. Group Contributions to Lens Potentials 

The understanding of environment-related biases and 
uncertainties in lensing constraints on the masses and 
shapes of galaxy dark matter halos, Hq, and substruc- 
ture requires detailed lens modeling (sec KZ04). That is 
beyond the scope of this paper and will be treated sepa- 
rately. For our purposes here, a simple way to quantify 
environmental effects is to determine the dimensionless 
convergence k and shear 7 that the group contributes to 
the lens potential. The convergence represents gravita- 
tional focusing created by additional mass at the posi- 
tion of the lens galaxy, while the shear represents tidal 
distortions created by having an inhomogeneous distri- 
bution of matter near the lens galaxy. Convergence can 
never be constrained using lens models alone because of 
the mass sheet degeneracy (|Gorenstein et al1ll988l: iSahal 
I2OOQ1) . As a result, it is often omitted, which can lead to 
significant biases in the model results (KZ04). Shear can- 
not be constrained in models of two-image lenses, which 
leads to enormous model uncertainties. While shear can 
be constrained in models of four-image lenses, one of the 
puzzling results is that shear is required in nearly all four- 
image lenses, and the required shear strengths cannot 
easily be explained by tradit i onal models of large-s cale 
structure l|Keeton et alJll997t IPalal fc Watsonll2004^ . 

Models of four-image lenses lead to the rule of thumb 
that a shear of 7 ~ 0.1 is common for groups, and 7 ~ 
0.3 for clusters. Generally, we expect k > 7 because of 
the way convergence and shear add when there is more 
than one perturber (see the Appendix). Many lensing 
conclusions scale as (1 — k) to some power (see KZ04 and 
the Appendix), which can help us estimate the biases. 
For example, if a lens has convergence k ~ 0.1, then lens 
models that omit the convergence will overestimate Hq 
by ^^10%. Hence, we consider convergences and shears 
larger than ~0.1 to be quite important, and values down 
to ~0.05 worth consideration. 

In this section, we quantify the effect of the group envi- 
ronments in our sample by determining the convergence 
and shear due to the group surrounding each lens galaxy. 
Because it is not clear how the mass is divided between 
the individual group members and a common group dark 
matter halo, we consider two extreme cases that bound 
the range of possibilities. In the "group halo limit," we 
assume that all of the mass is associated with a common 
group halo, and we estimate the amount of mass from 
the velocity dispersion. In the "group galaxies limit," we 
assume that all of the mass is bound to the individual 
member galaxies, and we calibrate the individual shears 
and convergences wi th recent weak lensing observations 
(jSheldon et al J 120051) . Comparing results from the two 
limits can indicate how much our conclusions depend on 
how the mass is distributed within the group. 

4.4.1. Group Halo Limit 

The group halo limit is appropriate for considering re- 
laxed, dynamically evolved groups in which the individ- 
ual dark matter halos of member galaxies may have been 
stripped by interactions, so that the dominant compo- 
nent is a common dark matter halo. In this case, the 



luminous galaxies just trace the underlying mass distri- 
bution, and the velocity dispersion of the galaxies is a 
measure of the mass of the group as a whole. If we model 
the group halo as a singular isothermal sphere (SIS), then 
we can take the velocity dispersion cr,. of a group and its 
centroid position (6, cf)) in polar coordinates centered on 
the lens galaxy, and compute the convergence Kgrp and 
shear 7grp using eq. (|A20I) in the Appendix. (Note that 
'^grp = 7grp in the SIS approximation.) We then use 
the uncertainties in the velocity dispersion and centroid 
position to determine the uncertainties in the conver- 
gence and shear. These derived uncertainties are highly 
non-Gaussian, so we compute them using Monte Carlo 
simulations. 

The results are presented in Table[5|* In three of the six 
groups (MG0751, PG1115, and B1422), the convergence 
and shear are significant. In these systems the group 
environment needs to be accounted for in lens models, 
and our measurements of the group properties make that 
feasible. In the other three groups (BRI0952, MG1654, 
and B2114), it appears from our current data that the 
convergence and shear are small. 

The case of B1422 illustrates how our identifications of 
additional group members affect estimates of the conver- 
gencc and shear. Based on six members, iKundic et alJ 
p997b) estimated that the group velocity dispersion is 
cTr — 550 km s~^, and that the centroid lies a projected 
distance of 6 ~ 14" from the lens galaxy. That led them 
to estimate a large "observed" shear of 7obs = 0.23, a 
value consistent with the shear 7mod ^ 0.2 required by 
lens models ( Hogg fc Blandford 1994; Kecton ct al. 199^ 
iDobler fc Keetoil2005l[l . Now with 16 members, we find a 
smaller velocity dispersion ar — 470 km and a larger 
offset b = 43", which reduce the nominal observed shear 
to a more moderate value of 7obs = 0.058. Accounting for 
the measurement uncertainties (which has not been done 
before), we find allowed ranges of 0.034 < 7obs < 0.095 
at Icr, and 0.020 < 7obs < 0.170 at 2a. It is clear that 
omitting the group members can bias the estimates of 
how the environment affects the lens potential. However, 
even with an extensive catalog (16 members), uncertain- 
ties are important when comparing observations of the 
group with inferences from lens models. The large dif- 
ference between the nominal observed shear 7obs — 0.058 
and the shear required by models 7mod ~ 0.2 might not 
be significant at more than 2a. To test consistency be- 
tween the observations and models, it is best to make 
new models with the environment constrained by our 
observations (including the uncertainties, and possible 
dumpiness within the group; see §4.1). 

4.4.2. Group Galaxies Limit 

In the group galaxies limit, we suppose that all the 
mass in the group is bound to the individual member 
galaxies. This is probably a better approximation for 
dynamically younger groups, which may still be in the 
process of collapse and whose galaxies still retain their 
halos. 

To calibrate the shear and convergence from each 
galaxy, we turn to observations that are perfectly suited 

* As noted in the table caption, we assume Zs = 2 for the un- 
known source redshift in B2114, but the particular value has little 
effect on our results. 
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TABLE 5 

Convergence and Shear in the Group Halo Limit 



Lens 


b Ur 


ft^grp — 7grp 


0^ 


Number of 




["] [km s-i] 




[dcg] 


Images 



PG1115 


23 


440 


0.089 


+0.0C5 
-0.046 


6 ±37 


4 


B1422 


43 


470 


0.058 


+0.037 
-0.024 


-65 ± 20 


4 


MG0751 


23 


320 


0.049 


+0.071 
-0.033 


-27 ±36 


4+R 


MG1654 


65 


200 


0.007 


+0.011 
-0.005 


-35 ± 30 


R 


BRI0952 


61 


170 


0.005 


+0.013 
-0.004 


-81 ±42 


2 


B2114'' 


49 


110 


0.003 


+0.014 
-0.002 


-12 ±52 


2 



Note. — The groups are sorted by decreasing values of the con- 
vergence and shear. The magnitudes of the convergence and shear 
are equal (Kgrp = 7grp) in the SIS approximation. Values larger 
than 0.05 are marked in boldface. The angle 9^ defines the direction 
of the shear (measured North through East). The convergence and 
shear errorbars are Itr uncertainties derived from the uncertaintines 
in the group centroid position and velocity dispersion (from Table 
Bl . Column 6 lists the image configuration for each lens: 2-image, 
4-image, or ring (R). 

''We assume a source redshift Zs = 2 for B2114, but the particular 
value has little effect on our results. In particular, assuming 2s = 3 
leads to the same numerical values for the Kgrp, 7grp, and B-f. 



to our needs: weak lensing. The advent of large surveys, 
such as the Sloan Digital Sky Survey, has made it possible 
to measure shear as a function of distance from an aver- 
age galaxy, on scales from 20 h^^ kpc to 7 h^^ Mpc (e.g., 
iSheldon et al.ll2005j) .^ The projected offsets of galaxies 
in our sample range over ~ 20-700 kpc, so we are 
workin g in precisely the regime studied by Sheldo n et al.l 
()2005tl . One possible concern with this calibration is 
that the Sheldon analysis provides limited information 
about how sh ear and co nvergence scale with luminosity 
(or mass). .Sheldon et al . (2005) measured the shear pro- 
file in three luminosity bins, but most of our galaxies fall 
into just one of the bins (—22 < Mi — 51og/i < —17). 
Thus, we expect that our present analysis characterizes 
the average properties of the group galaxies well, but it 
will be worthwhile to redo the analysis when future weak 
lensing data provide finer luminosity resolution. 

The quantity that is measured in weak lensing studies 
is shear, but the observed shear profile is consistent with 
a power law which means that there is a simple relation 
between the shear and convergence. The Appendix gives 
this relation, and also provides more details about how 
we use the weak lensing data to calibrate our analysis. 

Once we have computed the convergence and shear 
from each galaxy , we combine them in the manner of 
eqs. (jA12|) - (IA16|) in the Appendix. The net convergence 
and shear for each group are presented in Table We 
also list statistical uncertainties derived (using Monte 
Carlo simulations) from t he errorbars on the we ak lens- 
ing parameters quoted bv ISheldon et alJ l)2005j) . These 
are generally small and are probably less important than 
systematic effects due to incompleteness (see below). As 
in the group halo limit, in the galaxies limit we find that 
MG0751, PG1115, and B1422 aU have significant group 
contributions to the lens potential. The group contribu- 

^ ISheldon et alJ 1200,51^ quote comoving distances, but we have 
converted to angular diameter distances as those are more natural 
for our analysis. 



tions are fairly smaU in MG1654, BRI0952, and B2114. 
Comparing Tables |S1 and leads to a crucial point: our 
conclusions about which groups significantly affect the 
lens potentials do not depend on assumptions about how 
the mass is distributed within the groups. While the 
detailed lensing implications do depend on the mass dis- 
tribution — which is actually useful (see §4.4.3) — it is 
reassuring to see that our qualitative conclusions about 
which groups are important for lensing are robust. 

In carrying out this analysis we do not distinguish be- 
tween different morphological galaxy types. Given their 
higher mass-to-light ratios, elliptical (red) galaxies would 
produce more shear and convergence than spiral (blue) 
galaxies of the same luminosity. To test the effects of 
morphology on our results we consider the extreme as- 
sumption t hat all group galaxies are red and use the 
lSheldon~ al. (2005) shear profiles for red galaxies. In al- 
most all cases the total convergence and shear due to the 
group increase by ~70%. Thus, if most group members 
are red galaxies then our current analysis may actually 
underestimate the convergence and shear by as much as 
70%. 

As a simple sanity check, in Table |H1 we also compute 
the net convergence and shear assuming that each galaxy 
can be treated as a singular isothermal sphere (SIS) with 
velocity dispersion a — 100 kms^^. We expect that 
many of the galaxies have larger velocity dispersions, and 
given that convergence and shear scale as ct^, this sim- 
ple case should provide a conservative lower bound on 
the net convergence from the group. (The shear is more 
complicated, as discussed below.) Comparing the results 
assuming a = 100 km s~^ with those based on the weak 
lensing calibration confirms this expectation, and gen- 
erally suggests that the weak lensing calibration is rea- 
sonable. We should note that the convergence from the 
simple SIS model may not be a strict lower bound if ha- 
los are significantly truncated, but in practice most of the 
convergence (and shear) arise from galaxies close enough 
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TABLE 6 

Net Convergence and Shear in the Group Galaxies Limit 







[Sheldon] 




[a = 100 km s"^] 


Lens 


Ktot 


7tot 


[deg] 


Ktot 7tot 

[deg] 



64.2 ±3.2 0.040 0.031 66 

-35.61^3 5 0.064 0.019 -34 

-82.9 ± 1.8 0.040 0.029 -88 

21.4 ± 13.2 0.013 0.004 43 

-64.6 ±8.3 0.006 0.003 -67 

-15.9 ±3.3 0.008 0.007 -16 



Note. — The groups are sorted as in Table ISl Values larger than 0.05 are 
again marked in fjoldface. Columns 2-4 list the net convergence and shear 
when we calibrate the individual galaxies using the weak lensing observations 
by [Sheldon ct al. 12005). The errorbars represent la statistical uncertainties 
derived from the weak lensing errorbars quoted by Sheldon et al.; even more 
important may be systematic effects due to spectroscopic incompleteness (see 
text). For a simple comparison, Columns 5-7 list the results when we treat each 
galaxy as an isothermal sphere with velocity dispersion a = 100 km s~^. 

''Again, we assume a source redshift Zs = 2 for B2114, but the particular value 
has little effect on our results. 



PG1115 0.066 to ooe 0-028 003 

B1422 0.092 t^fol 0.010 t^H 

MG0751 0.115 too26 O-O^l ^0 034 

MG1654 0.027 toMi 0-006 to'oos 

BRI0952 0.013 to 002 0.007 15!'°°^ 

B2114^ 0.016 0.012 



to the lens that truncation -would have Httle effect. For 
example, if all halos had cut-off radii at 300 kpc, the SIS 
model convergence would drop from k = 0.040 to 0.036 in 
PG1115, and from 0.064 to 0.057 in B1422. Even with 
an extreme cut-off at 100 kpc, the convergences would 
still be 0.027 for PG1115 and 0.051 for B1422. (In aU 
cases, the shears are basically unchanged.) We expect 
that the conservative assumption of a small velocity dis- 
persion more than compensates for the omission of a cut- 
off radius, so that the SIS model results are indeed lower 
bounds. 

It is again interesting to consider how our efforts to in- 
crease the group membership have affected conclusions 
about the convergence and shear. The most instruc- 
tive case is PG1115. With our catalog of 13 members, 
we find a net convergence k = 0.066lQoJg and a net 
shear 7 = 0.028j:i^;^^j| at position angle 9^ = 64 ± 3. If 
iKundic eTa l. (1997a) had done the same analysis with 
their catalog of four members, they would have found 
K = 027g °°i 7 = 0.014lS]-ol]6, and 9^ = 69 ± 5. 
ijTonr-yl [l998l would have obtained similar results.) In 
other words, the previous catalogs missed at least half of 
the sources of shear and convergence. The problem was 
that they focused on the region within ^^30" of the lens, 
but (in the galaxies approach) the galaxies have extended 
dark matter halos, requiring the inclusion of group mem- 
bers out to the full virial radius of the group in order to 
fully characterize environmental effects in lens models. 

With this thought in mind, we must consider how spec- 
troscopic incompleteness may affect the results in Table 
El We cannot account for galaxies that we missed, but 
we can ask how our results would have differed had we 
omitted a few of the galaxies that we did actually in- 
clude. The effects of incompleteness on the net conver- 
gence are simple: since convergences from diff erent galax- 
ies combine in a simple scalar sum (see eq. IA12I in the 
Appendix), omitting galaxies causes us to underestimate 
the net convergence. Turning this around, we can say 
that our estimates of Ktot are strict lower bounds on the 



convergence from the group. To be more quantitative, 
in analyzing subsamples of our group catalogs we find 
that the net convergence scales roughly with the num- 
ber of group members; so if the true number of members 
is, say, 50% larger than what we have observed then we 
expect the true convergence to be ~50% larger than our 
estimate. For understanding biases in lens models, it is 
very valuable to have a lower bound on the total conver- 
gence, because that can be turned into lower bounds on 
the biases. 

The shear is more complicated, because multiple 
cont ribut ions sum as tensors rather than scalars (see 
eqs. EH and Em ; adding more contributions can ei- 
ther increase or decrease the net shear, and modify the 
position angle. The effects depend on the spatial dis- 
tribution of member galaxies. In PG1115 the most im- 
portant galaxies lie roughly in a line on the sky, which 
means that the direction of the net shear is robust against 
incompleteness, while the amplitude of the shear scales 
roughly with the number of group members. In contrast, 
in B1422 the galaxies are distributed more broadly, so in- 
completeness may change the shear direction by tens of 
degrees and the shear amplitude by tens of percent. (Of 
course, the shear from the galaxies in B1422 is small, so 
even large fractional uncertainties are not so important.) 
Finally, in JVIG0751 the environmental effects are domi- 
nated by the galaxy Gl that is massive and close to the 
lens. 

The bottom line is that incompleteness does not sig- 
nificantly affect our qualitative conclusions. Since our 
estimate of the net convergence is a lower bound, we 
know that our conclusion that groups are important on 
MG0751, PG1115, and B1422 is robust. While it is pos- 
sible that our conclusion that the groups are not so im- 
portant in MG1654, BR10952, and B2114 could change 
if we measure more galaxies and find that k rises, the 
observed correspondence between th group velocity dis- 
persion and richness (see §4.2) suggests that adding more 
members would not affect k significantly. Incompleteness 
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issues will need to be considered when making detailed 
quantitative comparisons between our environment ob- 
servations and lens models. 

4.4.3. Discussion 

It is important to understand the similarities and dif- 
ferences between the results from the "group halo" and 
"galaxies" limits. We have already noted that both limits 
lead to the same conclusions about which groups arc im- 
portant for lensing. All three high-dr groups (MG0751, 
PG1115, and B1422) produce a significant convergence 
in both approaches, so these groups cannot be ignored in 
lens models. At the same time, two of the low-Ur groups 
(BRI0952 and B2114) produce a small shear and conver- 
gence in both approaches, suggesting that these groups 
are not so important for lensing. While this latter conclu- 
sion may not seem exciting, it is actually quite valuable. 
Two-image lenses (including both BRI0952 and B2114) 
suffer from a strong degeneracy between ellipticity and 
shear, if both quantities are unknown. That degeneracy 
can now be broken by ruling out models with large shear. 
The situation is less clear for MG1654, because the group 
halo analysis implies negligible shear, while the galaxies 
analysis yields a small but non-negligible shear 7 = 0.03 
(and that could be an underestimate). 

When we turn to a more quantitative comparison of the 
two approaches, we notice some significant differences. 
The differences suggest that it may be possible to dis- 
tinguish between the group halo and galaxies limits, and 
thus to learn about the distribution of dark matter within 
the groups. For example, in B1422 the group halo analy- 
sis leads to a moderate convergence and shear, while the 
galaxies analysis leads to a larger convergence but a neg- 
ligible shear. Lens models require a large shear 7 ^ 0.2 
that is marginally consistent with the group halo ap- 
proach (given our uncertainties) but grossly inconsistent 
with the galaxies approach. In PG1115, as in B14 22, lens 
models require large external shear 7^0.1 (Jmpev et al.l 
IT998I) . The group halo limit result (7 = 0.089t[!!]|^) is 
within 1(7 of the model requirement while the galaxies 
limit gives a factor of three lower shear, inconsistent with 
predictions. The results in both B1422 and PG1115 — 
the two highest velocity dispersion groups — suggest that 
the mass is distributed in a common halo rather than be- 
ing attached to the individual group members (as might 
be expected if high velocity dispersion groups are more 
dynamically evolved). This hypothesis needs to be ex- 
amined more carefully with detailed lens models; rather 
than just comparing the shear required by lens models 
with that inferred from our observations, it is important 
to build models that explicitly incorporate the environ- 
ment (which may even consist of multiple subgroup halos; 
see §4.1) in both the group halo and galaxies limits and 
see whether either case can fit the lens data. Systems like 
B1422 and PG1115 may provide an exciting opportunity 
to determine the distribution of dark matter in a distant 
group. 

We also note that KZ04 used PG1115 as a fiducial 
example with w hich to asses envi r onme nt-related biases 
in lens models. iDalal fc WatsonI l)2004f) suggested that 
PG1115 is a very atypical lens environment, and that 
KZ04 therefore overestimated environmental effects. To 
the contrary, we find that PG1115 is nothing if not typi- 
cal: the group's kinematic properties and shear and con- 



vergence are consistent with at least half of the groups 
in our sample. 

Another interesting system is MG0751. Here, the 
galaxies analysis is dominated by the Gl galaxy, ly- 
ing just 6" from the lens. Even so, Lchar et al. (l99'3) 
showed that lens models including only the lens galaxy, 
Gl, and up to three other nearby galaxies cannot fit the 
lens data. It will be interesting to use new lens mod- 
els to test the hypothesis that both Gl and the common 
group halo contain significant mass, and to see whether 
we can constrain their relative masses. We must issue 
two warnings, however. First, the projected offset of Gl 
from the lens galaxy is just 20 kpc, which lie s at th e 
inner limit of the range studied bv lSheldon eTall (|200l|); 
thus, the reliability of the weak lensing calibration is not 
clear. Second, we show in §4.5.2 that MG0751 also has 
a significant shear from a group along the line of sight, 
which must be included in lens models along with the 
group at the lens redshift. 

In summary, our results above show that: (1) Group 
environments, whether the mass lies with individual 
member halos (galaxies limit) or in a common group halo 
(group halo limit), can contribute significantly (k,7 > 
0.05) to the lens potential. (2) If the members have ha- 
los (galaxies limit), they can have a big effect, perhaps 
even greater than that of a common group halo. (3) In 
the galaxies limit, correcting for incompleteness is only 
going to boost the convergence (but will move the shears 
in either direction). (4) In the highest velocity dispersion 
groups, the shears produced in the group halo limit are 
more consistent with the observationally required values, 
suggesting that we might be able to discriminate between 
the models of the mass distribution. 

Finally, it is remarkable that the three lenses with sig- 
nificant environmental effects include both quad lenses 
(PG1115 and B1422) and the one quad/ring (MG0751), 
whereas the double/ring (MG1654) and the two double 
lenses (BRI0952 and B2114) all have small convergence 
and shear. While this result is limited by small number 
statistics, it may suggest a correlation between image 
confi guration and environment. Conventional wisdom 
(e.g.. iRusin et alJl200l() holds that shear does not sig- 
nificantly affect the relative numbers of quads and dou- 
bles. However, KZ04 argue that treating environment 
properly (including terms beyond a simple shear) does 
change the quad/double ratio. Our new results provide 
empirical evidence that there is a connection. 

4.5. Lensing Effects of Line-of-Sight Structures 

Since lensing is a projected phenomenon, we must 
consider whether structures projected along the line 
of sight significantly affect strong lens systems. Dif- 
ferent theoretical approaches to studying the effects 
of interlopers on strong lensing in a ACD M universe 
have yielded contradictory result s (e.g., | Bar-Kana| 
1996; ' Kceton et al. 1997; Prcm adi fc Mart cl 200i 
l^ambsganss. Bode fc Ostrikeil 1200^) 7" so an empir- 
ical approach is clearly necessary. To date, there 
are only three lenses with confirmed line-of-sight 
groups (B0712-h47 2, B160 8+656 and M G 1131 - ^0456 ; 
[Fassnacht fc LubinI 120021: Tonrv & Ko chanekl 120001: 
fFassnacht et al. 2005) and sev eral other candidates 

( Faure et al 2004: .Morgan et al.i2005) . Here we present 

the first systematic survey for structures along the 
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line-of-sight to strong lens systems. The large redshift 
baseline and wide field of view of our spectroscopic 
survey make it ideally suited to address this issue. 

4.5.1. The Zone of Influence 

Our first task is to estimate the "zone of influence" for 
a perturber along the line of sight to each lens in our 
sample. Although it has not been done before, the cal- 
culation is straightforward using the formalism presented 
in the Appendix. Briefly, if we assume that a perturber 
can be modeled as an isothermal sphere with some given 
veloci ty dispersion, then we can use eqs. (|A6|) . IIA7p . and 
HA20|I to compute the effective convergence and shear 
(kbS and 7eff) as a function of the impact parameter b of 
the perturber relative to the lens, and the redshifts of the 
lens galaxy and perturber. (In the SIS approximation, 
KeS = 7cff for a single perturber.) 

Figure El shows contours of k^s — 7eff in the plane of b 
and Az = Zpcrt — ^icns, for perturbers with velocity dis- 
persions of 100, 300, or 500 km s~^. (Because we do not 
see rich clusters along the lines of sight to these lenses, 
this range of ar should span the observed range of struc- 
tures.) The grayscale is explained in the figure caption. 
We consider the zone of influence to be the region in 
which Koff,7cff > 0.05, i.e., the unshaded (white) region 
in the figure. The shape of this region is not sensitive to 
the fact that an SIS halo has an infinite extent; a cut-off 
halo radius of 300 kpc does not change the zone of infiu- 
ence. Also interesting are regions in which Keff and 7eff go 
negative (shaded black), which represents a breakdown 
of our formalism. This happens only when the intrinsic 
convergence of the perturber is k > 0.5, which means 
that the offset between the lens and perturber is small 
enough that the lens actually lies within the Einstein ra- 
dius of the perturber (see the Appendix). In this case, 
the "perturber" is no longer just a perturbation because 
its caustics interact with those of the main lens galaxy, 
and we would observe a strong lensing effect from the 
second mass as well. This breakdown does not affect our 
conclusions because we do not actually see any structures 
lying within this region; besides, any objects that lie so 
close to the line of sight to the lens would presumably be 
known from previous observations. 

As the velocity dispersion of the perturber increases, 
the zone of its influence grows dramatically oc a^. Con- 
sequently, a more massive perturber can produce a large 
shear and convergence even when offset from the lens; a 
perturber with ctj. 500 km (i-e., a rich group) can 
be offset by as much as 1' and still produce 7 ~ 0.05. 
Another striking feature of Figurel^lis the very wide red- 
shift baseline in front of and behind the lens over which a 
perturber can cause significant effects. This result shows 
that it is crucial to catalog not just mass structures in 
the immediate vicinity of the lens galaxy, but also else- 
where along the line of sight, in order to model the lenses 
accurately. 

4.5.2. Prominent Interloping Structures 

To quantify line-of-sight effects for the lenses in our 
sample, we consider the most prominent structures 
(those likely to be groups and clusters) identified from 
our redshift catalog. This approach is analogous to the 
group halo limit for groups around lens galaxies (§4.4.1), 



For figure see file f6_sm.jpg 

Fig. 6. — Contours of KcH ~ 'yaB produced by a perturber 
with a velocity dispersion of 100, 300 and 500 kms~^, lo- 
cated at a redshift < Zpat < 1 and having an impact pa- 
rameter b with respect to the lens, computed for seven of 
the lenses in our sample. (We exclude PMN2004 because 
its lens redshift is unknown.) The horizontal axis represents 
the redshift difference between the perturber and lens galaxy. 
The vertical axis represents the projected distance of the per- 
turber from the lens, out to 2'. Contours are drawn at k^b — 
7eff = 0.001,0.01,0.05,0.1 and 0.5, although not all of them 
are clearly visible in all panels. The grayscale is as follows; 
^^011,708 > 0.05 (white); 0.05 > Koff,7cff > 0.01 (light gray); 

0. 01 > Kcff, 7cft > 0.001 (medium gray) and fi:cff,7off < 0.001 
(dark gray). Important regions are 0.5 > Kcff,7cff > 0.05, 

1. e., the areas in white. At small impact parameter, the per- 
turbation approximation breaks down as we enter a regime in 
which the perturber itself causes strong lensing (see text); in 
our calculations, this leads to negative values of k^b and 7cfi 
(shaded in black). Notice the strong dependence on a: the 
zone of influence scales as a^. Massive perturbers can have 
large effects even when they lie far from the lens. Another 
striking feature is the very wide redshift baseline in front and 
behind the lens over which the perturber can cause a signifi- 
cant effect. 



because it accounts for the dark matter in massive bound 
structures. To identify potentially important line-of- 
sight structures, we show the redshift histograms again in 
Figured now shading only those galaxies that lie within 
1' of the lens. We have over-plotted the cur ve of the nor- 
malized shear strength Kcs/k = 7eff/7 fea. lAll)) to give 
an indication of how the convergence and shear vary with 
redshift. From the discussion in §4.5.1, we expect struc- 
tures with large velocity dispersions, small projected off- 
sets from the lens, and/or small redshift offsets from the 
lens galaxy to contribute most to the lens potential. To 
be conservative, we select only those peaks in Figured 
that: (1) have at least four members; (2) lie within Az 
such that Kcs/k and jcs/"f are > 0.5; and (3) have at 
least one member projected within 1' of the lens. For 
every peak, we set pessimistic 3<t velocity limits and use 
bi- weight estimators of location and scale to calculate the 
mean velocity Vi and the line-of-sight velocity dispersion 
(Jr.i- Based on the membership, we then calculate the 
projected spatial centroid of the structure, and its offset 
from the lens g alaxy bj an d posi tion angle (pi. Finally, 
we can use eqs. (|A6|I . (|A7|I . and ljA20|l to determine the 
effective shear and convergence. 

The results are presented in Table d We list the in- 
dividual effects of all prominent structures. We choose 
to include the groups at the lens redshifts (repeating re- 
sults from Table [SJ, so that we can compare local versus 
interlop ing structu res. We then sum all structures (us- 
ing eqs. IA12HA16)l to obtain the total convergence and 
shear for each lens (in the group halo limit). We empha- 
size that our results are conservative in the sense that we 
have not tallied all the line-of sight objects that might 
affect the lens models. We do not include galaxies out- 
side prominent peaks and plan to address this question in 
future work. We have not included all peaks in the veloc- 
ity histograms (Fig.di because most are undcrsampled 
or intrinsically poor, and thus have too few members for 
us to interpret them as likely groups and to compute a 
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Fig. 7. — Same as the left panels in Figure U] except that here we have shaded only those galaxies within l' of the lens. We have 
also overplotted the normalized shear strength and normalized convergence {k^s/k = 7ofi/7; see text) to guide the eye regarding general 
behavior of convergence and shear as the perturbing structure is moved away from the lens in redshift. 



meaningful velocity dispersion. In addition, the veloc- 
ity dispersions we compute are probably underestimates 
not only because of the narrow (conservative) choice of 
initial velocity ranges but also because many of the veloc- 
ity peaks are poorly sampled. Finally, our spectroscopic 
target selection prioritizes galaxies thought to lie at the 
lens redshift (see §3), and to some extent that limits our 
ability to identify interloping structures. 

The discussion of incompleteness in §4.4.2 applies here 
as well, with one minor change. As before, we note that 
convergence is a scalar and that the contribution from an 
overdensity is always positive, and conclude that iden- 
tifying additional prominent structures would only in- 
crease the total convergence Ktot- The change is that 
an uncertainty in the convergence zeropoint prevents us 
from declaring that we have obtained a strict lower bound 
on the total convergence from the line of sight. The zero- 
point uncertainty arises from the effects of underdensities 
such as voids. Strong lensing calculations conventionally 
assume that there are a few density peaks superposed on 
top of a smooth background at the mean density of the 
universe. Voids can then be thought of as regions where 
the density is negative (relative to the mean), which con- 



tribute nemiiwe convergence to the lens potential (e.g., 
ISeliaklHool . If there is a significant negative conver- 
gence, it effectively changes the convergence zeropoint: 
in summing the effects of overdense structures along the 
line of sight, we should start not from zero but from the 
appropriate negative convergence. 

We cannot measure the convergence zeropoint directly 
because that requires knowledge of the local density of 
matter (dark and luminous) at every point along the line 
of sight, which is impossible to obtain even with complete 
redshift surveys. Nevertheless, we can make a useful es- 
timate based on a simple model. We create Monte Carlo 
simulations of random lines of sight in a universe in which 
some of the mass is contained in halos while the rest is in 
a smooth background at a level below the mean density. 
In these simulations we are able to determine the total 
convergence Ktot from density fluctuations along the line 
of sight, as well as the contribution Kpeaks from promi- 
nent structures (^Cpcaks is analogous to what we compute 
above from our data). We can then interpret the differ- 
ence, = Ktot — Kpeaks, as thc Convergence zeropoint. 
(Details of the calcu lation, and further d iscussion of the 
results, are given bv lKeeton et"ani2nn5bl 'l 
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Fig. [7] — continued. 



Considering many random lines of sight, we find that 
Kzp has a mean of zero and a dispersion of < 0.02, for typ- 
ical lens and source redshifts and reasonable halo mass 
functions. Because the zeropoint uncertainty is small, 
we believe that our Ktot values above are probably lower 
limits on the total convergence. In other words, it is un- 
likely that voids contribute enough negative convergence 
to counter the positive k from prominent peaks in a given 
line of sight. 

The main result is that four of the eight lenses 
in our sample have significant interloping structures. 
MG0751, PG1115, and B1422 each have one structure, 
and HE2149 has three. At least one of those structures, 
along the line of sight to MG0751, has a significant con- 
tribution to the lens potential. This perturbing group lies 
at a redshift of z = 0.56, which places it between the lens 
galaxy and the source quasar, and has a velocity disper- 
sion Ur — 550 km s~^ based on the six members that we 
have identified. With a centroid that lies at a small pro- 
jected offset of 28" from the lens, the group contributes 
a convergence and shear kios — 7ios — 0.066. This is 
clearly an important contribution to the lens potential 
— in fact, it is slightly stronger than our estimate of the 
contribution from the group at the lens redshift (in the 



group halo limit). It represents one more piece in the 
interesting puzzle of fully understanding lensing in the 
MG0751 system. 

To our knowledge, our lens sample is not biased toward 
having significant line-of-sight effects. Our survey meth- 
ods are, if anything, somewhat biased against finding 
line-of-sight structures (as discussed above). Therefore, 
our discovery of a significant structure in 1 /8 lenses sug- 
gests that line-of-sight effects are important in at least 
~10% of all lenses. That estimate needs to be confirmed 
with a larger sample, but it does indicate that lensing 
effects from the line of sight deserve further attention. 

5. CONCLUSIONS 

We have presented the first results from our spectro- 
scopic survey of the environments of strong gravitational 
lenses. We have used multislit spectroscopy to measure 
the redshifts of 355 galaxies in the fields of eight strong 
gravitational lenses with lens galaxies at redshifts be- 
tween 0.25 and 0.50. After adding 16 redshifts from the 
literature, we have analyzed a total sample of 371 galax- 
ies with redshifts. 

The lens galaxy belongs to a poor group in six of the 
eight systems in our sample. We discover three new 
groups associated with the lens galaxy of BRI0952 (five 
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TABLE 7 

Convergence and Shear, Due to Prominent Line-of-Sight Structures 



Lens 


ID 


^pcrt 


b 

["] 


[km s^-"-] 


l^cB = 7eff 


0^ 
[deg] 


Ktot 


7tot 


9j,tot 
[deg] 


MG0751 
MG0751 
MG0751 


1 
2 

total 


0.35 
0.56 


23 
28 


320 
550 


0.049 
0.066 


-27 
-56 


0.116 


0.101 


-44 


BRI0952 
BRI0952 


1 

total 


0.42 


61 


170 


0.005 


-82 


0.005 


0.005 


-82 


PG1115 
PG1115 
PG1115 


1 
2 

total 


0.31 
0.49 


23 
44 


440 
300 


0.089 

0.010 


7 

-74 


0.099 


0.080 


5 


B1422 
B1422 
B1422 


1 
2 

total 


0.34 
0.28 


43 
80 


470 
400 


0.058 

0.020 


-65 
-28 


0.078 


0.066 


-57 


MG1654 
MG1654 


1 

total 


0.25 


65 


200 


0.007 


-35 


0.007 


0.007 


-35 


HE2149 
HE2149 
HE2149 
HE2149 


2 

3 

4 

total 


0.27 
0.45 
0.60 


51 
70 
56 


400 
180 
150 


0.017 
0.004 
0.003 


-62 
-80 
88 


0.024 


0.022 


-68 


B2114 
B2114 


1 

total 


0.31 


49 


110 


0.003 


-12 


0.003 


0.003 


-12 



Note. — "Prominent" structures are defined as having at least four members, at 
least one of which is projected within 1' of the lens galaxy and displaced by Az such 
that Keff/K = 7cff/7 > 0.5. Column 2 labels all the prominent structures along the line 
of sight to each lens; "0" refers to a group at the lens redshift. (HE2149 is the only lens 
with more than one interloping structure.) Columns 3-7 refer to individual structures, 
while Columns 8-10 give the final results after combining all the prominent structures 
along the line of sight to each lens (including a group at the lens redshift, if there is 
one). Values larger than 0.05 are marked in boldface. 



members), MG1654 (seven members), and B2114 (five 
members). We more than double the number of mem- 
bers for another three previously known groups around 
the lenses MG0751 (now 13 members), PG1115 (13 mem- 
bers), and B1422 (16 members). These six groups add 
to the still small number of all poor groups identified at 
intermediate redshifts. 

We determine the kinematics of the six groups, includ- 
ing their mean velocities, velocity dispersions, and pro- 
jected spatial centroids. For the newly discovered groups, 
we quantify these properties for the first time. For the 
other three groups, the increased membership allows us 
to make more robust estimates of the kinematic proper- 
ties of the groups than previously possible. The highest 
velocity dispersions we measure (320 to 470 km s~^ for 
MG0751, PG1115, and B1422) are consistent with those 
of nearby dynamically-evolved. X-ray luminous groups 
(MZ98), while the lower velocity dispersions (110 to 200 
km s~^) are more typical of dynamically younger groups 
at low redshift. In the two cases where a diffuse X- 
ray component ha s been measured (PG1115 and B1422; 
iGrant et al.ll2003) . the X-ray temperatures and our ve- 
locity dispersions are consistent with the local cr^-Tx re- 
lation (MZ98). 

To understand the evolution of groups and their galax- 
ies, it is important to determine the relation of the bright- 
est group galaxy (BGG) to the group potential. (In 
four of the six groups, MG0751, BRI0952, PG1115, and 
B1422, the lens galaxy is not the BGG.) We find that 
the BGG generally lies off the center of the group po- 
tential and occupies an orbit indistinguishable from the 
other group members. This result is surprising in com- 
parison with nearby. X-ray luminous groups, in which 



the BGG is always a giant elliptical galaxy occupying 
the center of the potential, with an orbit distinct from 
the other group members (ZM98). However, most of the 
effect we see comes from the three groups with lower ve- 
locity dispersions. In two (MG0751 and PG1115) of the 
three highest velocity dispersion groups, the BGGs lie 
within the errors of the group centroid, suggesting that 
at least these systems are comparable to dynamically- 
evolved poor groups in the local universe. 

We use our detailed observations of the groups to as- 
sess how environments affect gravitational lens models. 
A key ingredient is an accurate determination of any off- 
set between the lens galaxy and the group centroid on 
the sky. In MG0751, PG1115, B1422, and MG1654, the 
lens galaxy is offset spatially from the group centroid. 
Obtaining a larger sample (which is underway) to deter- 
mine the full distribution of lens vs. group offsets will 
be important for understanding how lens environments 
affect statistical quantities such as the quad/double ratio 
and lensing constraints on Ha (see KZ04). 

To quantify environmental contributions to lens poten- 
tials in more detail, we estimate the convergence (gravi- 
tational focusing) and shear (tidal distortions) from each 
group. We consider two different models of the group 
mass distribution that bound the extremes of dynami- 
cal states of groups. The members of young groups are 
likely to still have large dark matter halos and the group 
mass may be dominated by the dark matter halos of the 
individual member galaxies. In this approach, we cali- 
brate the shear and convergence from each galaxy based 
on observations of weak lensing in the Sloan Digital Sky 
Survey (Sheldon et a l. 2005). and then sum the contri- 
butions appropriately. As the group evolves, these halos 
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may be truncated via interactions, so the mass will be 
redistributed into a common group halo that may be the 
dominant group mass component. In this approach, we 
approximate that halo as an isothermal sphere and use 
the measured group centroid and velocity dispersion to 
compute the convergence and shear. 

At least three of the lenses in our sample (MG0751, 
PG1115, and B1422) have convergences and shears large 
enough (ac,7 > 0.05) to indicate that the environment 
plays a significant role in the lens potential. For these 
systems, our survey substantially improves the observa- 
tional constraints that will be needed to make detailed 
lens models that properly include environmental effects. 
Remarkably, the high shear and convergence values occur 
in the quad lens systems, while the environments of the 
double lenses are relatively weak. This result suggests 
that environment may affect the relative numbers of quad 
and double lenses, a topic much debated in the literature 
Csee IKing et al.l Il9 96: Kochanek 1996 bt IKeeton et all 
1997^ 'Rusin & Tegmark 200j ICohn fc Kochanekl 120041 
Keeton fc Zabludofi 2004). For the other lenses, the con- 
clusion that environment does not significantly affect the 
lens potential is also valuable: constraining previously 
unknown environmental terms to be near zero will still 
improve lens models. 

For the first time, we present a systematic assessment 
of whether structures along the line of sight to lens sys- 
tems are important for lensing. We show that interlop- 



ing structures can in principle affect lens models over a 
wide range of spatial and redshift offsets. Our pencil- 
beam survey is ideally suited to identifying such struc- 
tures if they are present. We find that at least four out of 
eight lenses have prominent linc-of-sight structures, i.e., 
groups whose spatial and redshift offsets place them in 
the "zone of influence" of the lens. MG0751, PG1115, 
and B1422 each have one substantial group along the 
line of sight, while HE2149 has three groups at differ- 
ent redshifts. Of these, the interloping group in MG0751 
has a significant effect on the lens potential. Our sur- 
vey is actually biased against interloping groups (and is 
not complete), so finding that at least one of eight lenses 
10%) is affected by projected structures is intriguing 
and worth further study. 
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APPENDIX 

FORMALISM FOR CONVERGENCE AND SHEAR 

IKeetoiJ l|2003f) presents a formalism for deriving a simple analytic estimates of the effective convergence Kcs and 
shear 7efi produced by a perturber somewhere along the line of sight to a lens. The same formalism can be used for 
all the situations considered in the main text (see §§4.4 and 4.5). The only assumption is that the convergence and 
shear are small, so that we can work at first order in both quantities. We will show below that there is a context in 
which this assumption breaks down, but with no effect on our conclusions. 
If the perturber were at the same redshift as the lens galaxy, it would produce a convergence and shear given by 
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where (^pcrt is the lens potential of the perturber. Here 7 is the shear strength, and 6^ the shear direction (which we 
measure North through East). If t he per turber lies at a different redshift Zpcrt 7^ zi, then the convergence and shear 
are modified to the effective values ()Keeton..2003.') 



Keff = 



7off = 



(1-/3) [n-P{K^-^^)\ 



(1 -/3k)2 - 
(1-/3)7 



where 



D{zi,Z2) 



D{Q,Zs) 



D{0,Z2) D{zi,zs) 
zi ==min(z;, Zpcrt) , 
Z2 = max(zi,Zpert) , 



(A6) 



(A7) 



(A8) 

(A9) 
(AlO) 



Lens Galaxy Environments 



23 



where D{zi,Z2) is the angular diameter distance between redshifts zi and Z2- Note that when Zpert = zi, we have 
[3 — Q and so we recover KcS — k, and 7eff — 7. Moving the perturber away from the lens redshift increases (3 and 
decreases KcS and 7eff. 

In certain circumstances, the effective convergence and shear can apparently go negative, which seems puzzling. 
The sign flip occurs only when the denominator in eqs. ljA6|l - ljA7|l goes negative. This can happen only when n and 
7 are sufficiently large — in particular, only when the line of sight passes through the strong lensing regime of the 
perturber. (For an isothermal sphere [below], this corresponds to k,7 > 0.5.) In this case, the lensing critical curves 
of the perturber would merge with those of the main lens galaxy, which would completely change the configuration 
of lensed images. In other words, our formalism breaks down when the "perturbation" is sufficiently strong, but if it 
were that strong it would (presumably) be known already. 

If K and 7 are small, then we can expand eqs. HA6|I - (|A7I) to first order write 

M^M^l.^. (All) 

K 7 

We can think of Kcs/k and 7off/7 as the "normalized" convergence and shear — the actual perturbation strength, 
normalized by the value that would apply if the perturber were at the same redshift as the lens galaxy. We use this 
quantity in the text as a simple way to characterize the redshift dependence of the perturbation strength. 

We now specify how to determine the net convergence and shear when multiple perturbers are present. This is 
relevant when we assess the effects of a group around a lens by considering the member galaxi es (§ 4. 4.2) , and also 
when we consider the effects of structures along the line of sight (§4.5.2). We see from eqs. (IA1|) - (|A3|) that the 
quantities k, 7c, and 7^ are linear in the perturber potential. The shear is a combination of 7^ and 7^ that actually 
corr esponds to a rank-2 tra celess tensor, or a headless vector (headless because it is invariant under rotation by 180°). 
fSee iSchneider et alJll993 for more discussion.) Thus, if we know the effective convergence KoS.i and shear 7cff,i, as 
well as the shear position angle O^^i (North through East), for a set of perturbers, then the proper way to compute 
the net effects is as follows: 

Ktot = ^ Kcff,,; . (A12) 

■i 

7c,tot = ^ 7eff,i cos 26*-^,^ , (A13) 

i 

7s,tot = ^ 7cff,i sin 26l^,i , (A14) 



7tot = y7?,tot +7^,tot , (A15) 
1 



''7, tot - 



tan ^(7s,tot/7c,tot) • (AI6) 



Note that KcS^i > for any real perturber, so the terms in the convergence sum all go in the same direction. By 
contrast, the cosine and sine factors mean that terms in the shear sums may add or cancel. Hence, incompleteness in 
our sample can only cause us to underestimate Ktot, but it may cause us to over- or underestimate 7tot- 

It is worthwhile to recall how convergence and shear affect lens models. Convergence is largely responsible for 
systematic biases in lens models (KZ04), through the mass sheet degeneracy ( Gorcnstcin et al. 1988: Saha 2000,1 . The 
biases can be thought of as simple rescalings of model parameters, such as: 

/3K(l-«:tot), (A17) 

/icx(l-Ktot), (A18) 

fix{l-Ktot)-\ (A19) 

where /3 is a mass parameter related to the lens velocity dispersion, h is the Hubble parameter, and fi is the image 
magnification (see KZ04 and references therein). Neglect of convergence is the main source of biases in lensing results 
for quad lenses. Double lenses, by contrast, are so under-constrained that poor knowledge of convergence and shear 
leads to lens models that are just plain wrong. In both cases, detailed observations of lens environments are necessary 
to derive the constraints necessary to make lens models reliable. 

While the formalism presented so far is fully general, it is valuable to discuss two particular perturber models. First, 
if we can approximate a perturber as an isothermal sphere then we can easily relate measurable quantities to k and 
7. Specifically, for an isothermal sphere with velocity dispersion a and impact parameter b, we have 

We can measure the line-of-sight velocity dispersion (Tr of the group and the offset b between the lens galaxy and group 
centroid, use them to determine k and 7, and finally fold in the redshift difference as above to determine and 7cff. 

The second specific model we consider is a power law density profile calibrated by weak lensing. In a large and 
detailed analysis of galaxy-galaxy weak lensing in the Sloan Digital Sky Survey, iSheldon et al.l l)2005(l present the 
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average shear as a function of radius for a sample of 127,001 lens galaxies. They find that the shear profile is consistent 
with a simple power law 

7(i?) = i?"" , (A21) 

■^crit 

where Scrit = {c^ Dos) I {^t^GDoiDis) is the critical surface density for lensing, which carries all the dependence on the 
lens an d source redsh ifts through the angular diameter distances Doi, -Dos, and Dis between the observer, lens, and 
source. iSheldon et al.l tabulate the values of a and A for galaxies in three luminosity bins.^ For any given galaxy that 
we observe, we convert from our measured J-b and apparent magnitude t o absolute magnitude in the SDSS z-band, 
using the color and iiT-corrections computed with lBruzual fc CharlotI ljl993D spectral synthesis models. (The luminosity 
bins used by Sheldon et al. are wide enough that small systematic uncertainties in the magnitudes do not shift galaxies 
between bins.) We then look up the values of A and a for each luminosity from Table 2 of Sheldon et al. We use the 
lens and source redshifts to compute Scrit and hence ^{R)- The last thing we need is the convergence. For a power 
law 7(i?) oc there is a very simple relation between shear and convergence: 

n(R) = ^— ^ 7(i?) . (A22) 
a 

Thus, it is straightforward to compute k and 7 for each galaxy. We can then factor in the redshift distance relative to 
the main lens galaxy as described above. 



® ISheldon et al.l 120051) actually tabulate power law parameters for the galaxy— mass correlation function rather than the shear directly, 
but their formalism makes it straightforward to convert back to A and a. At any rate, those parameters are more fundamental in terms of 
what they measure. One additional technical point is that Sheldon et al. quote lengths using comoving distances, but we convert those to 
angular diameter distances in our analysis. 
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